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ABSTRACT 
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The coupling of N — 1 superconformal field theories of type (4m, 2) to two- 
dimensional supergravity can be formulated non-perturbatively in terms of a dis- 
crete super-eigenvalue model proposed by Alvarez-Gaume, Itoyama, Manes and 
Zadra. We derive the superloop equations that describe, in the double scaling 
limit, the non-perturbative solution of this model. These equations are equivalent 
to the double scaled super- Virasoro constraints satisfied by the partition function. 
They are formulated in terms of a c = 1 theory, with a Z2-twisted scalar field and 
a Weyl-Majorana fermion in the Ramond sector. We have solved the superloop 
equations to all orders in the genus expansion and obtained the explicit expres- 
sions for the correlation functions of gravitationally dressed scaling operators in 
the Neveu-Schwarz and Ramond sector. In the double scaling limit, we obtain 
a formulation of the model in terms of a new supersymmetric extension of the 
KdV hierarchy. Integrating over the fermionic eigenvalues we obtain a representa- 
tion of the discrete supersymmetric free energy in terms of the free energy of the 
one-matrix model. 
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1. INTRODUCTION 



1.1. String Theory and 2D-Quantum Gravity 

String theory [1] is our most promising candidate to be the theory that unifies 
quantum field theory with general relativity. The classical motion of a string in 
.D-dimensional space-time is described by a surface parametrized by a function 
X^(a,r) with /i = 1,...,D, which depends on two coordinates, a and r. To 
quantize the theory we are interested in the calculation of the following functional 
integral: 



Z = J VXe~ s ^ x \ (1.1) 

where the integration is done over all possible surfaces (space-time trajectories of 
the string). Here S[X] is the action functional. The definition (1.1) of the partition 
function is a natural generalization of Feynman's path integral for point particles. 
For a string propagating freely in flat Minkowski space, the action is determined 
by the proper area swept by the string and is given by the expression 



S[X] = T J dadr^\detg ab \ = T J dadr\J {X ■ X') 2 - X 2 X' 2 (1.2) 



where 



Here the proportionality constant T is the string tension, (^°,^ 1 ) = (r, a) and the 
Minkowski metric is rj^ = diag(l, — 1, . . . , — 1). 
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This is known as the Nambu-Goto action. However, it is awkward to work with 
this action because of the square root and of the fact that it is highly nonlinear. 
A more convenient form that is classically equivalent to the previous one can be 
obtained if we introduce an auxiliary field g a b, that is the metric tensor of the 
string world-sheet 

S[g,X] = ~J dadr^-gg a % v d a X^d b X\ (1.4) 

where g = det^. Polyakov used this action to formulate string theory in terms 
of the path integral 



Z = J VXVge~ S[9 ' X] . (1.5) 

The action (1.4) is conformal invariant (or local Weyl invariant), which means that 
its form does not change under a rescaling of the metric: 



9ab - e^g ab . (1.6) 

This invariance only holds in two dimensions since in this case \f^g is proportional 
to while g ab is proportional to e"^. The energy-momentum tensor is obtained 
from the variation of the action (1.4) with respect to the metric 



5S 
^5g~ a 



Tab = -7^=7-5, (1.7) 



and is given by the expression: 
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T ab = d a Xd h X - l -g ah g cd d c Xd d X = 0. 



(1.8) 



This energy-momentum tensor vanishes identically because of the equation of mo- 
tion of the field g. It is therefore also traceless. This tracelessness reflects the 
fact that the action is invariant under conformal transformations of the metric. 
The actions (1.2) and (1.4) are both invariant under arbitrary reparametrizations 
(a, t) — > (f a (cr, r), f T (cr, r)) of the surface, under which 

9ab - d a f%f l g kl . (1.9) 
The equation of motion for the field X^ is 

~^=fa {^gg ab d h X^ = AX» = 0. (1.10) 
For a free closed bosonic string we impose a periodic boundary condition 



X^(r )( j + 7r) =X fl (r,a). (1.11) 

The simplest topology of the world-sheet is a cylinder (Fig. 1) parametrized 
by (a, r), with < a < n and — oo < r < oo. 

We can use the reparametrization (diffeomorphism) invariance to fix = e^r} a bi 
where r\ TT = 1, i] aa = — 1, r\ Ta = and e v is an unknown conformal factor. In this 
(conformal) gauge the action for the field X fi (z,z) is 

S = ~\J d 2 ad a X^d a X^ (1.12) 
and the equation of motion takes the form 
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Fig. 1: Propagation of a free closed string. 



d 2 d 2 . 

' X» = 0. (1.13) 



dr 2 da 2 

The conformal transformation z = e w , where 



w — -2i(r + a) and w = -2i(r-a), (1-14) 

maps the cylinder to the complex punctured plane (D* = C — {0}. The equation of 
motion for X^(z, z) is then 



d z d s X"(z,z) = 0. (1.15) 

The most general solution compatible with the boundary conditions (1.11) is given 
by the following expression: 



We now introduce two different approaches to quantize the free bosonic string. 
The first one is the old covariant approach, based on a description in terms of the 
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coordinates only, and with restrictions to the physical Fock space in terms of the 
Virasoro conditions. This quantization procedure ensures that we have no ghosts, 
i.e. states with negative norm in the spectrum. It is similar to the Gupta-Bleuler 
quantization of the electromagnetic field. 

The energy-momentum conservation: 

d z T- z - z = d- z T zz = (1.17) 
implies that we can make a Laurent expansion around z = z = 0: 

T(z)=T zz (z)=J2LnZ~ n - 2 , (1.18) 

and the same for the anti-holomorphic part. Here L n and L n are the Virasoro 
generators that satisfy the Virasoro algebra: 

{L m ,L n } = -i(m -n)L m+n , (1.19) 

where { , } is the canonical Poisson bracket of the classical theory. After quanti- 
zation the Fourier coefficients of X fl (z, z) become operators that satisfy the com- 
mutation relations 

[<,<] = -m<Wn^- (1-20) 

The classical algebra (1.19) becomes then the central extension of the Virasoro 
algebra 
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[L m , L n ] = {m- n)L m+n + —m(m 2 - l)5 n+m , 



(1.21) 



where D is the anomaly term that follows from (1.20). 

Previously we saw that the energy-momentum tensor vanishes classically. If 
we would like to keep this property also in the quantum theory then the physical 
states should be annihilated by all the modes of this field. However it turns out 
that this constraint is too strong and not compatible with the anomaly D of the 
Virasoro algebra (1.21). Therefore, the best we can demand is that physical states 
are annihilated by the positive-frequency part of the energy-momentum tensor: 

L n |phys) = L n |phys) for n > 

(1.22) 

(L - &)|phys> = (Z - 6)|phys) = 0. 

Here b is a constant that comes from the ambiguity in the normal ordering of Lq. 
It can be shown that the spectrum is ghost-free provided that 6=1 and D = 26, 
which is the so-called critical bosonic string [1]. Remarkably, analog Virasoro 
constraints will also play an important role in the formulation of two-dimensional 
gravity in terms of the one-matrix model, as we will see later. 

In the modern covariant quantization of string theory we consider the path 
integral (1.5). A two-dimensional closed oriented manifold is topologically com- 
pletely characterized by the genus or number of handles h of the surface, so that 
we can write (1.5) as follows: 



where A is a loop-counting parameter. For each topology we have to sum over all 
possible embeddings (the VX integral) and over all possible metrics in that fixed 
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topology modulo diffeomorphisms. This is why the volume of the diffeomorphism 
group has to be divided out. In general, the action appearing in the above ex- 
pression describes a D-dimensional bosonic string. We will consider the case of 
a two-dimensional bosonic string or more precisely of a minimal conformal field 
theory coupled to gravity. 




Fig. 2: Topological expansion of the partition function. 

To actually calculate the integral J Vg, we have to distinguish between criti- 
cal and sub-critical strings. While for critical strings the integration simplifies to a 
finite-dimensional integration, in the case of sub-critical strings we are left with the 
hard problem to know how to quantize Liouville theory. In two space-time dimen- 
sions this 'continuum formulation' of quantum gravity has been first developed by 
Polyakov [2] and by Knizhnik, Polyakov and Zamolodchikov [3] who formulated the 
quantization in the light cone gauge and obtained the gravitational dimensions of 
gravitationally dressed primary fields for spherical topologies. Using the conformal 
gauge, David [4] and Distler and Kawai [5] generalized their results to Riemann 
surfaces of higher genus. In general, we are interested in the calculation of string 
amplitudes with external string states, so that we have to insert vertex operators 
in the functional integral and we would like to calculate the correlation functions 
of these operators. Several methods have been suggested to compute correlation 
functions in Liouville theory coupled to minimal conformal models. These tech- 
niques allowed the computation of one-, two and three-point functions [6,7] and in 
special cases iV-point functions [8,9] for spherical topologies. 
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More recently Mukhi and Vafa [10,11] have shown that for the special case of 
the D — 1 model coupled to gravity it might be possible to obtain the amplitudes of 
tachyon vertex operators on surfaces of arbitrary genus directly in the continuum 
approach. This can be done using the equivalence to a topological field theory 
formulated in terms of a Kazama-Suzuki model or a topological Landau-Ginzburg 
theory [12]. 

Another possibility to evaluate the sum over random surfaces in two-dimensionall 
quantum gravity is to discretize the problem and to replace the integral over the 
metric in (1.23) by a sum over random triangulated surfaces of arbitrary topol- 
ogy. This approach reproduces (1.23) in the continuum limit. Surprisingly it has 
been found in 1989 that such a sum can be exactly evaluated using matrix model 
techniques [13]. 

The simplest matrix model is the one-matrix model, where the field is an 
Hermitian N x N matrix. Using the planar loop equations, Kazakov [14] has 
carried out the analysis of the critical regimes of these models. It has been realized 
that they describe the coupling of non-unitary minimal series of conformal field 
theories (CFTs) of type (2m — 1,2) coupled to two-dimensional quantum gravity 
[15]. 

With the introduction of the double scaling limit [13], it has been possible 
to evaluate these models non-perturbatively, i.e. on surfaces of arbitrary genus. 
It turns out that in this case there is a close relation to integrable hierarchies 
[16,17], which is manifest through the appearance of the Korteweg-de Vries hi- 
erarchy (KdV), which describes correlation functions of scaling operators. The 
partition function is related to the tau-function of the KdV hierarchy. It has been 
shown [18,19] that the non-perturbative loop equation that fixes the solution of the 
one-matrix model in the double scaling limit can be obtained using the KdV struc- 
ture and the string equation as initial condition. This string equation describes 
the form of the renormalized cosmological constant in terms of scaling operators 
at the k th multicritical point. The loop equation is equivalent to an infinite set of 
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constraints satisfied by the partition function. They form a Virasoro algebra and 
are therefore called Virasoro constraints. They can be found directly in the matrix 
model representation of the partition function [20,21] or can be elegantly derived 
in the Witten-Kontsevich theory [22]. These constraints also play an important 
role in the generalization of the bosonic theories to N = 1 supersymmetric CFTs 
coupled to two-dimensional supergravity, as we will later see. 

1.2. SUPERSTRINGS 

In the last section we have introduced the simplest string theory, i.e. the bosonic 
string. However, it has been realized that all phenomenological interesting string 
theories should incorporate supersymmetry. Therefore we have to look at more 
general string theories, the superstring theories [1,23], which include fermionic 
degrees of freedom and supersymmetry relates the bosonic space-time coordinates 
X fi (a, r) to the fermionic partners \I/ M (<7, r). String theories of this type have many 
advantages since undesirable features of the bosonic string can be eliminated in the 
supersymmetric case as, for example, the presence of the tachyon. The free closed 
superstring (without boundaries) is described by the action* 

S = -— [ dadr (d a X»d a X v - m e a d a ^ v ) r}^. (1.24) 
2tt J 

Here g a represent the Dirac matrices satisfying the two-dimensional Clifford alge- 
bra: 

{g a ,g b } = 2 V ab , (1.25) 

where r/ ab , is the flat Minkowski metric. The matrices g a can be represented as 
follows: 

* The string tension T is related to the string coupling constant a' by T = (2na')~ 1 . In the 
following we choose the convention a' = 1/2. 
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"J and (1.26) 



The world-sheet spinor has two components in this basis 




(1.27) 



Since g a are purely imaginary, the Dirac operator ig a d a is real and it makes sense 
to choose the components of the world-sheet spinor real. Such a two-component 
real spinor is known as a Majorana spinor. The symbol ^> means g°. 

The action (1.24) is invariant under the global transformations 



5x v = ^ 

(1.28) 

= -ig a d a X^e, 

where e is a constant anticommuting Majorana spinor. Since these transformations 
mix bosonic and fermionic coordinates they are called supersymmetry transforma- 
tions. If e is not constant, the variation of the action (1.24) is 



5S = ^ J dadr(d a e)J a . (1.29) 
Here J a is the conserved supercurrent 



J a = \Q b Qa^d b X^ (1.30) 
that satisfies d a J a = because of the classical equations of motion 
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d a d a X fl = and g a d a V = 0. 



(1.31) 



Introducing the coordinates (z,z) as in (1.14) and the notation = ^ and 
^ + = ^ the gauged fixed action takes the form 

S = -7^J d 2 z {d- z Xd z X - M^c^ - t^ttj . (1.32) 
The energy-momentum tensor and the supercurrent for the free superstring are: 



(1.33) 

1 

2 



Tf(z) = Jz{z) = —WdzXn 



and similarly for the antiholomorphic components. 

Supersymmetry can be nicely formulated in the superspace language. The 
world-sheet coordinates z, z are supplemented by a pair of anticommuting variables 
9, 9 so that we have the supercoordinates: 

z = (z,9), z = (z,9). (1.34) 

Taking into account the fermionic character of 9 and 9, any function /(z,z) can 
be written in the following way : 

/(z, z) = f (z, z) + 5)0 + f 2 (z, z)9 + / 4 (z, z)99. (1.35) 

The supersymmetric derivative is defined by D = d/d9 + 9d/dz and the integration 
rules are: 
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J do = o, J dee = i and J dedeee = 1. (i.36) 

Introducing the superfield 

X"(z, z) = X^, z) + 9^(z, z) + 0^(2, 5) + OOG^z, z), (1.37) 
where G tJ '(z,z) is an auxiliary field, we can write the action (1.24) as follows: 

S = -— [ d 2 zDX^DX^ (1.38) 
2ir J 

where the infinitesimal volume element in superspace is d 2 z = d 2 zdedB. The super 
energy-momentum tensor T(z) has two components: 

T(z) =T F (z) +eT(z) (1.39) 
and after quantization we obtain the following OPEs for T(z) and T F {z): 



T{Z)T{W) = Xz—^f + J^y2 + + 



rW3>W = |^ + p&0 + ... (1.40) 

2 (z — w) z (z — u>) 



T F (z)T F (w) - 



(z — w) 3 2 (z — w) 
here .D is the range of the index ji. 
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An important point when dealing with fermionic strings is the boundary con- 
dition. A closed bosonic string has periodic boundary conditions X^(r,a + tt) = 
X^(t,<j). For the superstring the same holds for X^, while for the world-sheet 
fermions we can have periodic or antiperiodic boundary conditions. The first 
one is called the Ramond (R) boundary condition while the second one is the 
Neveu-Schwarz (NS) boundary condition: 

#**(t, a + tt) = # m (t, a) (R) 

(1.41) 

tt"(r, a + n) = -^(r, a) (NS). 

Since ^(z) is a (1/2,0) field, the transformation properties under z — > e 2m z in 
the two sectors are 

q?»( e 2ni z) = (R) 

(1.42) 

^(e 2ni z) = V»(z) (NS). 
The fields ^^(z) can be expanded in both sectors in Fourier modes, as follows: 

___ ! { n e Z for (R) 

»<•<,) =5">~. * for (NS) (L43) 

The Neveu-Schwarz states therefore describe space-time bosons, while the Ramond 
states describe space-time fermions. The field X^(z, z) has the same expansion as 
in the bosonic string. 

If we are considering fermions on a torus, there exist four possibilities to com- 
bine the boundary conditions around the two non-contractible cycles 

(R,R) (NS,NS) (NS,R) and (R,NS). (1.44) 
More generally, for a Riemann surface of genus h there are 2 2h possible choices of 



13 



boundary conditions, since there are four possibilities for each handle. Each of these 
possible assignments is called the spin structure of S/j. The boundary conditions 
for the holomorphic and anti-holomorphic fields \l/ and \I/ have to be chosen in 
such a way that modular invariance of the partition function is guaranteed. In 
the partition functions of superconformal minimal models the holomorphic and 
antiholomorphic spin structures are identified [24]. 

To quantize the free closed superstring in the old covariant approach we can 
proceed in analogy to the bosonic theory: we impose the (anti-)commutation rela- 
tions 

(1.45) 

The Virasoro constraints are now replaced by super- Virasoro constraints. The 
energy-momentum tensor has an expansion of the form: 



T(z) = J2L n z~ n - 2 (1.46) 

where L n are the Virasoro generators. The supercurrent has the same boundary 
condition as the fermionic fields. Therefore it has the following Laurent expansion 



1 _ 3 rGZ for (R) 

T F (z) = -Y G r z- r -^ with <^ . (1.47) 

W 2^ r }rGZ+ for (NS) 1 ' 



where G r are the generators of the superconformal transformations. The (anti- 
commutation relations follow from the OPE (1.40): 
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[L m , L n \ = (m - n)L m+n + -m(m 2 - l)£ n+ m,0 

[L m , G r ] = (^— - rj G m+r (1.48) 

{G r , G s } = 2L r+s + - ^r 2 - 5 r+s fi, 

where c = D. In the Neveu-Schwarz sector the five generators Li, Lq, G_ 1 / 2 
and G4/2 form a closed subalgebra. 

The constraint equations on physical states have the form of super- Virasoro 
conditions. In order to derive them we have to generalize the action (1.24) by 
incorporating a zweibein e%{r,a) 

K(3 = ele b pr) ab , (1.49) 

and an auxiliary field Xa( r ? cr ) (gravitino), where the indices satisfy a,b = 1, . . . , D 
and a, (3 — 0, 1. The complete action is then formulated in terms of these fields 
[25,26] 



S = -^-J drdaVhrj^ {h^d a X^d p X v - ie^g a d a ^ 

- i J drdaVhr ] ^e P h XaQ h Q a ^dpX v (1.50) 



i- J dTdaVh V ^^e^ b Xag b Q' 



a xp- 



The energy-momentum tensor and the supercurrent can be obtained from the above 
action by using 
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(1.51) 



The equations of motion for the the zweibein and the gravitino yield the constraint 
equations: 

T aP = and J Q = 0. (1.52) 

The generalization of the conformal gauge to the supersymmetric case is the so- 
called superconformal gauge, where: 

e% = e<P5% and X a = Qai>- (1-53) 

In the classical theory we can still use a Weyl rescaling and super- Weyl transforma- 
tion to gauge away the remaining metric and gravitino degrees of freedom, leaving 
only e a a = 5% and Xa = 0. Taking into account this gauge condition, the desired 
constraint equations are 

T aP = daX^dpX^ - ^nx^x^ + U>» Qa d^^ + l -^ Qp d^^ = o 

(1.54) 

Ja = l^Qa^dpX^ = 0. 

These constraint equations are incorporated into the quantum theory by requir- 
ing that the positive modes of the energy-momentum tensor and the supercurrent 
annihilate the physical states 
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G r |phys) = for r > 



L n |phys) = for n > (1.55) 
(L -6)|phys> = 0. 

If the G1/2 and G3/2 constraints are satisfied, the rest of them follow from the 
above algebra. The constant b takes different values depending on whether we are 
in the Ramond or Neveu-Schwarz sector. It can be shown [1,23] that b = is a 
consistent choice for the Ramond sector. For the critical dimension D — 10, a 
decoupling of negative norm states takes place, if b = 1/2 for the Neveu-Schwarz 
sector. 

In the case of the bosonic theory we have already mentioned that the one- 
matrix model provides us with a powerful non-perturbative description. These non- 
perturbative effects are very important in the context of superstring theory. They 
are responsible, for example, for supersymmetry breaking. We are therefore faced 
with the question whether it is possible to find similar (non-perturbative) powerful 
methods as the one-matrix model to evaluate theories with N — 1 supersymmetry 
coupled to world-sheet supergravity. 

A first approach in this direction has been taken by Alvarez-Gaume and Manes 
in ref. [27]. Here a naive replacement of ordinary matrices by supermatrices 
has been proposed to simulate the sum over supergeometries. Taking (N\M)- 
supermatrices one is able to simulate iV-bosonic degrees of freedom and M fermionic| 
ones. However the result of this approach is disappointing because these models 
are completely equivalent to bosonic models based on (N — M|0) matrices. 

From the continuum point of view, Di Francesco et al. [28] have analysed the 
connection between certain unitary superdiscrete series of minimal models coupled 
to 2D-supergravity and the Manin-Radul hierarchy [29] . This is one of the known 
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super symmetric extensions of the ordinary KP-hierarchy, as we will see in detail 
later. However, they came to the conclusion that odd flows are inconsistent with 
the string equation, so that the relation to superintegrable hierarchies remains 
obscure. 

The most promising approach to the problem of finding a model that describes 
N — 1 superconformal matter coupled to 2D-supergravity has been found by 
Alvarez-Gaume et al. [30]. The guiding principle in ref. [30] was a set of super- 
Virasoro constraints satisfied by the partition function, which code the superloop 
equations. Since the Virasoro constraints [18,19] played a prominent role in the 
Witten-Kontsevich theory [22] of the one-matrix model, it is reasonable to expect 
that a similar set of super- Virasoro constraints should capture some important 
features of the supermoduli space of super- Riemann surfaces. The model obtained 
in this way is an 'eigenvalue' model. It is formulated in terms of a collection of 
iV-even and iV-odd eigenvalues (A^, 9i). 

It has been shown that in the continuum limit this model reproduces the string 
susceptibility and the scaling dimensions of (4m, 2)-minimal superconformal mod- 
els coupled to 2D-supergravity [30]. Correlation functions of gravitationally dressed 
scaling operators in the Neveu-Schwarz sector agree with those computed in the 
continuum super-Liouville theory in the planar approximation [8,31]. A compari- 
son between the discrete and continuum approaches can be found in ref. [32,33]. 
Very recently E. Abdalla and A. Zadra [34] have established a precise dictionary 
between the operators of the discrete super-eigenvalue model and those of the 
super-Liouville theory using similar methods as those developed by Moore, Seiberg 
and Staudacher for the bosonic model [35]. 

An evaluation of the planar model for arbitrary superpotentials has been car- 
ried out by Alvarez-Gaume et al. [36,37,38]. The consideration of arbitrary super- 
potentials is important in order to derive continuum super- Virasoro constraints. 
They are equivalent to the double scaled superloop equations. The solution of 
these equations up to genus two shows that the purely bosonic part of the model 
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(fermionic coupling constants set to zero) agrees with the one-matrix model. 

The Virasoro constraints in the double scaling limit turn out to be described 
in terms of a c = 1 theory with a ^-twisted scalar field and a Weyl-Majorana 
fermion in the Ramond sector*. The derivation of these constraints is the first 
step on the way to finding the integrable hierarchy of super- differential equations 
that describes the flows in the double scaling limit. 

The solution of the super- Virasoro constraints in the double scaling limit or 
equivalently of the superloop equations for arbitrary genus has been found in 
ref. [40]. The supersymmetric partition function can be expressed in a simple 
way through the partition function of the one-matrix model if the fermionic de- 
grees of freedom are integrated out. This property appears in the discrete matrix 
integral even before carrying out the continuum limit. This is important for the 
geometric interpretation of the model. The correlation functions of scaling oper- 
ators in the double scaling limit allow the study of the connection to the known 
superintegrable hierarchies. 

1.3. Outline of the Thesis 

In this thesis we will solve the N — 1 super-eigenvalue model proposed by Alvarez- 
Gaume et al. [30] to describe the coupling of superconformal matter of type (4m, 2) 
to 2D-supergravity. First we will solve the model in the genus-zero approximation 
using an arbitrary bosonic part of the superpotential and derive the expression for 
the scaling operators in the Ramond and Neveu-Schwarz sector of the theory. We 
compute the expressions for the scaling exponents and show that we get agreement 
with the continuum super-Liouville approach [41,42]. 

For arbitrary potentials a careful treatment of the continuum limit is in order. 
This will be done first for the case of the one-matrix model, which we feel has 
not been done in the literature explicitly enough. We are going to prove that the 

★ These constraints have also been derived in [39] . 
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super- Virasoro constraints in the continuum are described by a c" = 1 theory with 
a 22-twisted scalar and a Weyl-Majorana fermion in the Ramond sector. 

After deriving the corresponding continuum (non-planar) loop equations we 
obtain their solution up to genus two and show that the bosonic part of the model 
(fermionic coupling constants set to zero) agrees with the one- matrix model result. 

We then derive the non-perturbative solution of the model in the discrete and 
continuum theory. The partition function of the supersymmetric model turns out 
to have a simple representation in terms of the one-matrix model if the fermionic 
degrees of freedom are integrated out. This relation is also present in the continuum 
theory, as can be seen from the expressions for the scaling operators and the string 
equation. We will then analyse the superintegrable hierarchy behind this model. 

The content of the subsequent chapters is the following: 

Chapter 2: We introduce the (super)-Liouville approach to 2D-(super)-gravity. 
We compute the critical exponents and the string susceptibility in order to compare 
them with the discrete approach through the matrix model. 

Chapter 3: We explain in some detail the one- matrix model. In section 3.1 we 
show that the large- N expansion of the matrix model can be used to carry out the 
sum over topologies which we would like to carry out in string theory. In section 
3.2 we derive the discrete Virasoro constraints and the equivalent loop equations. 
The one-cut solution of the loop equations for arbitrary potentials and genus zero is 
explained in section 3.3. In section 3.4 we explain how the KdV hierarchy appears 
in the double scaling limit. Related to this hierarchy are the double scaled Virasoro 
constraints or equivalently the loop equation in the double scaling limit [18,19]. 

Chapter 4'- In section 4.1 we present the one- matrix model in terms of eigen- 
values from the point of view that can be easily generalized to the supersymmetric 
case. We present the model obtained in this way in ref. [30]. We derive the form of 
the superloop equations, for arbitrary genus in the discrete theory, that are equiva- 
lent to the discrete super- Virasoro constraints. In section 4.2 we present the planar 
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solution of the model for arbitrary superpotentials. We derive the expressions for 
the bosonic and fermionic loop operators in the planar geometry. 

Chapter 5: In section 5.1 we derive the continuum Virasoro constraints and 
show that the bosonic and fermionic loop operators become a Z2-twisted scalar 
bosonic field and a Weyl-Majorana fermion in the Ramond sector. We obtain the 
expressions for the scaling operators and show that the critical exponents and the 
string susceptibility agree with those of (4m, 2) SCFTs coupled to supergravity. In 
section 5.2 we solve the continuum superloop equations up to genus two and show 
that the bosonic part of the model is represented by the one-matrix model. In 
section 5.3 we obtain the non-perturbative solution of the supersymmetric model 
directly from the discrete eigenvalue integral, which has a simple relation to the 
one-matrix model. A similar relation to the one-matrix model appears in the dou- 
ble scaling limit. We derive the non-perturbative expressions for the correlation 
functions and analyse the properties of the hierarchy behind the model in section 
5.4. In section 5.5 we review the most important properties of the known super- 
symmetric extensions of KdV and compare them with the hierarchy that describes 
our model. 

Chapter 6: We present our conclusions. 
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2. PATH INTEGRAL APPROACH 
AND (SUPER-)LIOUVILLE THEORY 



In order to gain some basic ideas about the results we should expect in the discrete 
approach to 2D gravity and its supersymmetric generalization, we start with a short 
review about the path integral approach to 2D (super)gravity. Some good articles 
about the Liouville approach to 2D-quantum gravity are for example [32,43,44]. 

2.1. Bosonic String and Liouville Theory 

We begin with the bosonic theory. The partition function of a general conformal 
system with central charge c = D described by the action Sm is the following 
integral: 

Here we have divided out the volume of the diffeomorphism group as before. For 
a free bosonic string moving in a flat background the action has the form*^ 

S M = ^j d^g^daX^dbX^. (2.2) 

The cosmological term appearing in (2.1) is absent in the classical theory but in the 
quantum theory it is necessary for renormalizability [2] ; /iq is the bare cosmological 
constant. Geometrically A = J d 2 ^^fg represents the area of the surface. 

The integration measures in (2.1) must be defined more precisely. The T>X 
measure can be defined using the normalization of the Gaussian functional integral 

-k Here we have carried out an analytic continuation to the Euclidean signature in space-time 

iX° — > A" and on the world-sheet it — > r, so that g a b gets positive definite, 
f In the following we change the conventions for a': we will set a' = 2. 
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J V g 5Xe-W 5x ^ = 1 with 
The Vg measure can be defined similarly 



11**11; = J d\^-g5X5X. (2.3) 



J V g Sge-^ = 1 with ll^llj = J d\^(g ac g bd + 2g ab g cd )5g ab 5g cd , 

(2.4) 

where Sg represents a metric fluctuation at some point g a \, in the space of metrics 
of a genus-/;, surface. These measures are invariant under the diffeomorphism group 
(the group of reparametrizations) but not under conformal transformations of the 
metric g^ — > t^g a b- It turns out that 



V e , g X = V g X exp (J^$l{<P, gfj , (2.5) 



where 



S L (<p, g) = J ^-g ah d aV d bV + R V + fi c e^J (2.6) 

is known as the Liouville action. Here /i c is a (ultraviolet divergent) non-universal 
term. 

The space of metrics on a compact topological surface £ modulo diffeomor- 
phism and Weyl transformations is a finite dimensional compact space Aih, known 
as moduli space. If for each point r e Aih, we choose a representative metric 'gat,, 
then the orbits generated by the diffeomorphism and Weyl groups acting on g ab 
generate the full space of metrics on S. Thus given the slice g(r), any metric can 
be represented in the form 
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f*9 = e«g{r) 



(2.7) 



where /* represents the action of the diffeomorphism. Since the integrand of (2.1) 
is diffeomorphism-invariant, the result of the integral would be infinite, unless we 
gauge-fix the symmetry. This can be done by introducing a pair of anticommuting 
Faddeev-Popov ghosts (b, c, b, c), as in ordinary gauge theories: 



S gh (b, c, b, c,g) = J d 2 £^g (b zz V- z c z + 6 ff V z c*) . (2.8) 

The integration over the metric g can be written as an integration over the moduli 
[dr], an integration over the conformal factor ip and an integration over the ghosts 



Vg(gh) = VgbVgCVgbVgC. (2.9) 

The partition function that we would like to calculate can then be expressed as 



Z = J [dr]Vg V Vg{gh)V g X exp {-S M -S gh -^-J d\^j . (2.10) 

We can now choose a metric slice g = e^g. Under this transformation we already 
know that VX transforms as (2.5), while a similar result can be deduced for the 
ghosts 



V^gh) = V g (gh)exp \-—S L (<p,g)j . (2.11) 

The gauge fixing of the measure V g ip is complicated since the norm implicitly 
depends on the metric. In [4,5] the assumption was made that the Jacobian of the 
transformation 
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V eV ^pV eV ~{gh)V eV ~X = J(<p, g)V dV V~(gh)V d X (2.12) 
takes the form of a local Liouville action, so that the partition function becomes 

Z = J [dr]V^V d (gh)V d Xe- S (2.13) 

with 



S = S M (X,g)+S gh (b,c,b,c,g) + J d 2 ^ (a^c^c^ + bRip + /^) . (2.14) 



Here \x is the renormalized cosmological constant. The constants a, b and c are 
fixed by requiring invariance under 



9ab -»• 9abe a ^ 

(2.15) 

¥>(0->v(0-*(0- 



This invariance comes from the fact that Z is only a function of e^g = g. The 
result is 



~ 25 -D „ b . 

b=— , a=-. 2.16 

48tt ' 2 K ' 



After a rescaling </? — > a/12/(25 — -D) </?, which gives a canonical kinetic term, the 
energy-momentum tensor corresponding to this action has the form: 
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T = -i(^) 2 + |d 2 y, with g= V /^ y ^. (2.17) 

It has a central charge ci = 1 + 3Q 2 - We can now determine c in (2.13). Since 
(p has been rescaled, we will write this term of the action as e 7V3 and determine 7 
from the requirement that e 7< ^ behaves as a (1,1) conformal field. From the OPE 
with the energy-momentum tensor we obtain the conformal weight: 

A(e^) = A(e^) = -±7(7 - Q) = 1. (2.18) 

Therefore 



7 = -L (v / 25~^D - v/T^d) • (2.19) 



12 

Here the sign between the square roots has been chosen in such a way that the 
classical limit is recovered as D — > —00. From the Liouville part of the action 
(2.15) it becomes clear that h = 7 2 , so that the classical limit should coincide with 
7^0. From the above formulas it becomes clear that the domain D < 1 is the 
one where the Liouville theory can be most easily interpreted. Summarizing, the 
total action is 



S = S M (X, g) + S gh (b, c, b,c,g) + ^J (> b <9 aV 9<9 &V 9 + QRp + /ie^) , 

(2.20) 

with Q and 7 given by (2.17) and (2.19) respectively. 

An important quantity, that we now have to consider, is the partition function 
for fixed area. We can rewrite (2.10) for fixed genus h as an integral over A 
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l h = J dAZ h (A) exp (-£^) , (2.21) 



where 



Z h (A) = J V V VXe- s 5 (^J d 2 ^-A^j . (2.22) 

The integration over the moduli and the ghost determinant are included in T>X . 
The partition function for fixed area (2.22) has a scaling behavior as a function 
of the area that is universal. It is characterized by a critical exponent, called the 
string susceptibility T s t r (h) 



Z h (A) ~ A r ^- 3 . (2.23) 

Due to this behavior, it turns out that in some cases the complete partition function 
(2.21) is ill defined. This exponent can be determined from (2.22) using simple 
scaling arguments. Shifting ip — > ip + p/j, where p is a constant, and setting 
A = exp (p) we can deduce: 

Z h {A)=A-^ l - h ^h- 1 Z h {\). (2.24) 
Comparing with (2.23) we obtain 



T str (h) = 2-Q(l-h). (2.25) 

7 



For genus zero this gives 
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rV(O) = -L (d - 1 - v /(25- J D)(l- J D)) 



(2.26) 



We now consider a minimal CFT coupled to 2D quantum gravity. These models 
have a finite number of primary fields and are described by a central charge 



D = c PyQ = 1 



pq 



(2.27) 



where (p, q) are coprime integers with q < p. In general, the spinless primary fields 
of the matter theory exp(ipoX), with conformal dimension Aq: 



A = A£ s = — 



1 

4pq 



(rq — sp) 2 — (p — q) 2 



(2.28) 



where 1 < r < p and 1 < s < q, get dressed by the Liouville field after the coupling 
to gravity 



$ = e ipoA-(O e MO. (2.29) 

The value of (3 is obtained by demanding that (2.29) is invariant under conformal 
transformations, i.e. that $ is a (1, 1) field: 

Aq-±(P-Q)P = 1. (2.30) 
From (2.30) we obtain the relation 



, = ^EE£±^£±^ = 1 ±s j91 i l +2Ao . (2 . 31) 



Here the branch (3 < Q/2 reproduces the semiclassical limit, where (3 — > as 
D — > — oo. It has been argued in [43,44] that only these operators can exist. They 
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are usually called operators on the right branch. Nevertheless, operators on the 
wrong branch, that violate the above Seiberg bound have an interesting application 
in the context of Witten's two-dimensional black hole [45]. 

For a minimal model the value of T str (0) obtained from (2.27) is 



T st r(0) = * (2.32) 

If (p,q) — (m + l,m) where m — 3,4, . . ., i.e. if we consider the unitary models, 
we obtain T str {Q) = —1/m. However, for general (p,q) models the above value 
of T str (0) does not agree with the expressions obtained from the matrix models, 
where 



Tstr(0) = -. (2.33) 

p + q — 1 

It has been suggested in ref. [46] that this discrepancy can be solved if we measure 
the scaling of the partition function with respect to the operator of lowest dimension 
of the CFT. For the unitary minimal models this operator is the identity, while in 
the non-unitary case it has a negative dimension 



_ l-(p-g) 2 , u 

4pq 



and we obtain 



= (2^35) 

Pmin p + q-l 

The expression for the string susceptibility takes then the form (2.33) (for spherical 
topologies) if we replace Q/j in (2.25) by Q/(3 m i n . For the simplest non-unitary 
minimal models of type (2m — 1,2) this exponent has the value T str (0) = —1/m, 
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so that it coincides with the expression previously found for the unitary minimal 
models coupled to gravity. 

Apart from the partition function, we would like to compute amplitudes for a 
surface to pass a given set of points. For a surface with pinned points {X{\ this 
amplitude is given by 



G(X h ■■■,X N ) = (j[J ^(&)5 (X t - X(£i)) d%^, 



(2.36) 



where the average is taken w.r.t. the functional integral (2.1). Passing to the 
momentum representation this is equivalent to 



G(pi, ...,p N ) = (j[f \fme iPiX{ii) d 2 ^. (2.37) 

These are the correlation functions of (gravitationally dressed) scaling operators 
(2.29). As the partition function, these correlation functions have a universal 
scaling behavior as a function of the area. The scaling behavior of the one point 
function for fixed area 



MA) = ZjA) J V * VXe ~ S5 (/ fty/de™ - A ) J d 2 ^Vde ipoX e^, (2.38) 
defines a critical exponent called the gravitational scaling dimension A: 



F$(A) ~ A^ A . (2.39) 

With similar scaling arguments as the ones used to determine 7, it is possible to 
show that the one-point function has a scaling behavior as a function of the area: 
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F^(A) = A^F^l). 



(2.40) 



This gives A = 1 — p/j and therefore 



V25 - D - VT^D v 1 

Integrating over the area the one-point functions of scaling operators, scale as 
a function of the renormalized cosmological constant 



j dAZ h (A)F*(A) exp (-£-) ~ /,- r -+ A+1 , (2.42) 
o 

and the partition function behaves as 



Zh ~ ^ Vstr - (2.43) 

These representations are useful in order to compare with the results from matrix 
models. 

To quantize the theory we are interested in the calculation of iV-point correla- 
tion functions of scaling operators. This can be done in a flat background [9,7,6,8] 
or in a non-trivial background such as the 2D black hole [47] . We are not going to 
explain these results here. 
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2.2. Non-Critical Superstring 



In the following we consider the N — 1 SCFTs with central charge c = D [48], de- 
scribed by the action Sjn, coupled to two-dimensional supergravity. The partition 
function for a fixed topology of the base manifold is [42] 



Z h = J VE^VXe~ SM ^ E \ (2.44) 

spin 
structures 

Here we have used the superspace formulation introduced in [49,50,51], where 
means the super-zweibein. In this formulation the action of the N — 1 free massless 
superstring takes the form 



S M = -L / d 2 zEDXDX. (2.45) 

87T J 

The field E is the superdeterminant of the super-zweibein E^ and D is the co- 
variant super derivative. 

In the supersymmetric action the cosmological term of the form 



/i J d 2 zE, (2.46) 

is not needed for renormalizability, in contrast to the bosonic case. The only 
reason to include this term in the action would be that it is compatible with all its 
symmetries. 

The construction of the super-Liouville action is completely equivalent to the 
bosonic theory and has been carried out in [41,42]. First we have to define a 
reference frame E: 
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E = e^ L E. 



(2.47) 



After gauge fixing the action in the superconformal gauge it is assumed that Jaco- 
bian of the transformation of the measures takes the form of a local super-Liouville 
action: 

S S lI$l, E] = ^J d 2 zE (^-D a <$> L D a <$> L + QY^L^j , (2.48) 

where 



$ L (z, z) = ip(z, z) + 6^(z, z) + (hp(z, z) + 00F(z, z) (2.49) 

is the Liouville superfield, and ip and are the fields that we introduced in (1.53); 
F is an additional auxiliary field and Y is the supercurvature [49,50,51]. In the 
component language the above action is 

from which we obtain the energy-momentum tensor 

Tl(z) = —{dipdip - Qd\ - ^). (2.51) 

It has a central charge 

c L = 1 + 3Q 2 + i = ^ (1 + 2Q 2 ) . (2.52) 

From the gauge-fixing procedure we get a Faddeev-Popov superdeterminant 
that can be represented in terms of two ghost superfields : 
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B(z) = P(z) + 6b(z) and C(z) = c(z) + B 1 {z), (2.53) 



where (6, c) are fermionic ghost with spin (2, —1) and a central charge c\, c = —26, 
while (/?, 7) are the bosonic supersymmetry ghosts with a spin (3/2, —1/2) and a 
central charge c^ 7 = 11. 

The original partition function (2.44) is invariant under a simultaneous rescal- 
ing of the frame E and a shift in <&£, 

£ _^ ^(M)^ $ L _, $ L _ I a ( Z) z). (2.54) 

7 

This condition can be used to determine Q and 7. The requirement that this 
symmetry is preserved in the quantum theory is equivalent to the vanishing of the 
total central charge 

ctot = c + c gh + c L = 0. (2.55) 

The conformal anomaly of the matter and ghost system is c + c g h — c — 26 + 11. 
We therefore obtain 

Q = (2-56) 

The constant 7 from (2.47) can be fixed demanding that exp (7<&l) transforms like 
a density, i.e. like a super conformal field of dimension (1/2, 1/2). Therefore 

A(e^) = ~ 7 (7 " Q) = ^ (2-57) 
This equation has the solution 
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7 = I - \ = -j= (V9^d- VT=B) . (2.58) 

In the above expression we have taken the branch that reproduces the correct value 
of 7 in the semiclassical limit. We observe that the theory is well defined provided 
that c < 1. At c — 1, the theory undergoes a phase transition and the local ansatz 
for the Jacobian is no longer valid [42]. 

In order to compute the critical exponents we have to rewrite the partition 
function as an integral over L, that is the characteristic length scale of a super 
random surface: 



oo 

Z h = J dLZ h {L) (2.59) 
o 



with 



Z h (L) = J V1LV$ L 5 (^J d 2 zE - I?J e - s M-S gh -S S L_ 



(2.60) 



spin 
structures 



Here the dimension of L is that of the length because the Grassmann integration 
in L = J d 2 zE lowers the dimension by one. 

The string susceptibility is defined as: 



Z h (L) ~ L r ^( ft M. (2.61) 
Like for the bosonic string, it can be computed from simple scaling arguments: 
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F str (h) = 2 - (1 - h)Q = 2 + (c - 9 - v / (9-c)(l-c)) . (2.62) 

To determine the gravitationally renormalized scaling dimensions we have to 
distinguish between operators in the Neveu-Schwarz and in the Ramond sector. In 
(1.45) we saw that there exist 2 2h different spin structures for the matter theory, 
corresponding to the different boundary conditions of the fermions around the non- 
contractible cycles. Similarly, there exist 2 2h spin structures for the supergravity 
theory. World-sheet supersymmetry requires that each spin structure of the matter 
system be coupled to the same spin structure of the supergravity system [52] . The 
spin structure for the holomorphic and antiholomorphic fields are identified. 

The operators in the Neveu-Schwarz sector are similar to the ones of the bosonic 
theory and can be written in the form: 



^ NS = J d 2 zEexp(ikX + (3<5> L ) . (2.63) 

Since these operators are to be integrated they have to be (1/2, 1/2) fields. The 
on-shell condition is then 



Ao-§(/?-Q) = i. (2-64) 

We use the notation Ao for the scaling dimensions of the primary fields of the 
matter theory, exp(ik~K.). The above equation has the solutions 



/5=|±i v / Q 2 -4 + 8A . (2.65) 



We have to keep the operators satisfying (3 < Q/2. 
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The renormalized scaling dimension can be computed from the scaling of the 
two-point function (for fixed superlength) as a function of L, in analogy to the 
bosonic theory 



NS , P VT^- VI-C+I6A0 



A INS = 1 - - = 2 1 — — — -. (2.66) 

7 VT^c- 



The operators in the Ramond sector are more complicated, since they have spin 
[42,52]. The renormalized scaling dimensions in this case are given by 



aR y/l-c- V-c+ 16A 

A R = - — v — u . 2.67 



The minimal N — 1 SCFT are characterized by a pair of integers (p, q) with 
q < p and they have a central charge 



1 2(p-g) 



(2.68) 
pq 



These theories have a finite number of primary fields that have a conformal dimen- 
sion 



^. (^(.-^ izfcg. ( , 69) 



with 1 < r < p — 1 and 1 < s < q — 1. 

Unitary SCFTs correspond to (p, g) = (m + 2, m). The simplest SCFTs are 
those of type (p, g) = (4m, 2) with m — 1, 2, They have a central charge 
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These are the models that we are going to consider in the next chapters. The case 
c = corresponds to pure supergravity and it is the only unitary model of the 
above series. 

From (2.69) we obtain the conformal dimension for operators in the Neveu- 
Schwarz sector: 

A** 1 ' 1 = l{l ~ 2m + l) with < = 0,l,...,2m-l; (2.70) 
Am 

and in the Ramond sector 



(i + 1 _ rof (ro _ 1/2)2 i with 

Am 16 

(2.71) 

The string susceptibility can be obtained from (2.62). We have to take into 
account that if the theory is non-unitary it is natural to assume that the cosmo- 
logical constant couples to the lowest dimension operator [52]. Therefore, we have 
to replace Q/7 by Q/f3 m i n in (2.62), and obtain 



and Pmin = - — . (2.72) 



Pmin p + q-2 
Pmin is the dressing of the primary field with the lowest weight 



a mm 

A - — . (2.73) 

For the series of type (4m, 2) we obtain the same value of the string susceptibility 
as in the bosonic theory r str (0) = —1/m. 
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One of the advantages of the continuum formulation of 2D- super gravity is that 
most of the calculations done on the sphere can be carried out with methods similar 
to those of the bosonic theory. Correlation functions of scaling operators can be 
computed using an analogous Coulomb gas formalism. This has been done in 
[31,32] and the results can be taken as a guiding principle to find a more powerful 
discrete formulation valid for higher genus that reproduces the continuum results 
for spherical topologies. We are not going to enter in more detail into the continuum 
formulation of 2D-(super)gravity. Details can be found, for example, in the review 
[32]. 
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3. THE HERMITIAN ONE-MATRIX MODEL 



In chapter 1 we saw that in string theory we are interested in the evaluation of a 
sum over Riemann surfaces of all possible topologies. This can be done, using the 
so-called double scaling limit, in the context of matrix models as was discovered 
in 1989 [13]. The simplest case of a matrix model is the one-matrix model, where 
the field is an Hermitian N x N matrix. The parameter N will appear as a genus- 
counting parameter in the partition function and the double lines that appear in 
the Feynman diagrams will be useful to characterize the orientability of the surface, 
as we will later see. 

In this chapter we will derive the discrete Virasoro constraints and the equiv- 
alent discrete loop equation that determines the solution of the model. We will 
analyse carefully the case of arbitrary potentials and prove that in the continuum 
limit a set of two decoupled string equations appears. We will furthermore present 
the one-matrix model in the double scaling limit and show the appearance of the 
KdV hierarchy. We introduce some basic notions of completely integrable systems, 
such as the Lax-pair formalism for the KdV hierarchy and the generalization to 
the KP hierarchy through the formalism of pseudo-differential operators. 

We will concentrate in this chapter on those aspects of the one-matrix model 
that we consider relevant for the supersymmetric generalization. A more complete 
analysis can be found in some excellent review articles [32,53,54,55,56]. 

3.1. Topological Expansion as a Large- V Expansion of a Matrix 
Model 



The partition function of the one-matrix model is defined as the following matrix 
integral [57,58,59]: 



Here $ is an Hermitian N x N matrix and Ag is a free parameter related to the 




(3.1) 
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cosmo logical constant, as we will see later. The potential is defined as: 



V(*) = Y + Y,9k* k - (3-2) 

fc>3 

The integration measure on Hermitian matrices is: 

N 

d N ^ = \[d^ ii J] d(Re$y)d(Im$y). (3.3) 

i=l l<i<j<N 

This partition function can be evaluated perturbatively using Feynman rules. Ex- 
panding the exponentials, expression (3.1) is reduced to the calculation of Gaussian 
integrals: 



Z B = I d N $exp 



TV 



2A B 



« 2 )n(| o $(-3" P (-)-)- <-) 



The product runs over the different types of vertices. The propagator is the ele- 
mentary integral: 



tw-'^-f^r 81 -^. (, 5 , 

J^ 2 $exp(-2£tr$2) N **' 

and can be represented by a double line (see Fig. 3). 



i < i 

J * 3 



Fig. 3: The elementary contraction for hermitian matrices. 
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This double line notation was first introduced by 't Hooft to represent a gluon 
propagator in QCD. Since $ is Hermitian we have no twisted propagators and 
therefore the surface is orientable. More complicated graphs can be evaluated 
using Wick rules. Diagrammatically a vertex of type $ 4 takes the form shown in 
Fig. 4. 



1 J 



AT 



I 

n 



u n 



Fig. 4: The quartic vertex representing tr$ 4 



To make the above clear, we consider a particular example, where only is 
different from zero. The partition function has the expansion of the form 



N ^2 N ^4 

— tr$ 2 _ _, 4tr$ 4 

(3.6) 



N ( N \ 2 

= 1 " x-24(tr$ 4 > + ( -%4 J (tr$ 4 tr$ 4 ) + . . . . 



The simplest contribution, where n^ — l has the following form 
(tr$ 4 ) = (<^-><<I>^> + ($y$?*>(^S> + (^«^><^r*> 



A 2 A 2 
= {Sikbjjhih + SuS jk 5ji5 k i + 5u5ji5ji5 kk ) = (N 3 + JV + N 3 ) , 

(3.7) 

that diagrammatically can be represented as shown in Fig. 5. 
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Fig. 5: These are the three contributions to (tr$ 4 ) of (3.7). 

In order compare with the Liouville theory, we are only interested in the con- 
nected graphs of the matrix model. The generating functional for the connected 
diagrams is the free energy 

F = \ogZ B . (3.8) 

This is the quantity that has to be compared with the partition function of Liouville 
theory. 

In general, a term of order g p p in the expansion (3.4) defines a closed Feynman 
diagram with n p vertices of type g p . A graph with (713, . . . , n p ) vertices of type 
(<73, ... , g p ) will give a contribution: 

(g^ . . . g?)N'- p + v \ P B - V where V = £ n p , (3.9) 

to the free energy, since each vertex contributes a factor Ng p /A B , each of the P 
propagators a factor Ag/iV and each of the / closed loops (index loops) contributes 
a factor N, due to the corresponding index summation. Since p-lines are incident 
on the vertex of type p, the total number of propagators P of a closed graph is 
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Fig. 6: The solid lines represent a Feynman diagram of a 
$ 3 interaction, while the dual simplex is built out of triangles 
and represented by broken lines. 

2P = ^2,pn p . To see how a genus expansion of the free energy can be obtained we 
have to introduce the dual simplex associated to each of these diagrams. 

Every edge of the dual simplex corresponds to a propagator of the Feynman 
diagram, every p-gonal surface to a p- vertex and each vertex to a face (Fig. 6). We 
can therefore associate the total number of edges £, faces T and vertices V of the 
dual simplex to the quantities P, V and / of the Feynman diagram. 

For the dual simplex we know from the Gauss-Bonnet theorem that 

X = 2-2h = F-£ + V (3.10) 

holds, where X is the Euler-characteristic and h is the number of handles of the 
surface. 

In this way the matrix integral (3.4) automatically gives a random simplicial 
representation of surfaces. The large N expansion of (3.1) can be interpreted a 
topological expansion, as we will now see. 

With the previous notation (3.9) becomes 
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\\g n /N 2 - 2h ^ nM2 ^\ (3.11) 

p>3 



so that we arrive at the following expansion for F: 



OO 

F = J2N^ h F h with F h = Y,T^& np(P/2 ~ 1) U9?- (3-12) 

h=0 Sh V 

F/j is the free energy for fixed genus h. The sum in F^ runs over all graphs of 
the theory, which may have a symmetry, so that the order of the symmetry group 
Gh has to be divided out. From (3.12) we see that only the spherical topologies 
(h = 0) contribute in the large-iV limit. However, we are interested in summing 
over all topologies and it is at this point where the notion of double scaling limit 
becomes relevant. 

To make contact with the continuum Liouville approach, in particular to the 
partition functions (2.10), (2.21) and (2.22) we have to introduce the notions of 
the area and the cosmological constant in this context. We can imagine the surface 
built out of equilateral triangles with area 1/2. It is easy to see that the total area 
of the surface is given by 

p 

The free energy for fixed genus can then be written as a function of the free energy 
for fixed area 



F h = Yl A B F h(A) With = £' jgj^ II ( 3 - 14 ) 

where the prime in the sum means that the area is held fixed. 
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F = ^N 2 ~ 2h F h . (3.15) 

For later convenience we define 

A B = e-^°, (3.16) 
where no is the bare cosmological constant on the lattice. Therefore, 



F h = J2Fh(A)e-^ A . (3.17) 
A 

It turns out that the above sum on the lattice is ill defined [60]. This is because 
F(A), that is the number of different surfaces of a given area A, grows too fast as 
A — > oo: 



F h (A) ~ A T ^~ 3 exp (pL c A) , (3.18) 

where fi c is a cutoff dependent coefficient. We will later see that the universal 
critical exponent T s t r (h) coincides with the one introduced in the continuum theory 
(2.39) for the (2m —1,2) minimal models coupled to gravity. 

In the vicinity of the continuum limit the sum over A can be approximated by 
an integral: 

oo 

F h ~ J dAe-^ A F h (A) ~ (fi - ^ C ) 2 - T ^ h \ (3.19) 
o 

This describes the behavior of the free energy for fixed genus as a function of the 
cosmological constant, which is important in order to describe the double scaling 
limit. 
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We are able to keep contributions from Riemann surfaces of arbitrary genus 
in the sum (3.12), if the limits N oc and /j,q — > /j, c are correlated. This can be 
done by keeping the 'renormalized string coupling constant' n fixed as N — > oo 

k- 1 = n{hq - Aic ) 1 - r -(°)/ 2 . (3.20) 

This form of taking the continuum limit is known as the double scaling limit [13]. 
The full free energy has an expansion in powers of k : 

oo 

F = ]T K 2h ~ 2 F h . (3.21) 
h=0 

To define the continuum limit we have to introduce a constant a with the 
dimension of a length, that is the lattice cut-off. This means that the length of 
the basic triangulations are equal to a and no longer constant. The 'renormalized 
cosmo logical constant' t is defined according to: 

Aio - Ai c = a 2 t, (3.22) 

where t is kept fixed by tuning (/iq — /i c ) as a — > 0. Taking \i c = we obtain 

A B = 1 - a 2 t. (3.23) 

For the one-matrix model the value of the string susceptibility is T str = —1/m, 
so that we effectively have to keep 

K = a- 2 - 1 ' 171 ^ (3.24) 

fixed, as we will now see. 
* Here Fh has been rescaled. 
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3.2. Virasoro Constraints and Loop Equations in the Discrete The- 
ory 

The Virasoro constraints satisfied by the partition function were first discovered 
in the continuum theory [18,19]. These constraints are closely related to the in- 
tegrability and the KdV structure of the one-matrix model in the scaling limit. 
However, it has been realized that a similar connection to integrable hierarchies 
appears in the discrete theory. Here the integrable hierarchies are the Toda and 
the Volterra hierarchies [61,62]. The discrete Virasoro constraints can be deduced 
directly from the matrix integral, as we will now see. 

The infinite number of constraints satisfied by the partition function (3.1) 
follow by using the invariance of Z under a change of variables [20,21] 

$' = $ + e $ n+1 for n > -1, (3.25) 

where e is a real parameter. Here we have taken the analytical transformations 
into account, which are the ones we are interested in, as we will see in a moment. 
These transformations are classical symmetries of the potential that are equivalent 
to the equation of motion of <E> 

oo 

Y^kg k $> k - l = Q. (3.26) 

For the quantum theory we have to take into account the change of the measure. 
Under this transformation we have = $ fc + ek& n+k and this implies for the 
measure and the potential to first order in e^: 



f For n = — 1 the measure is invariant. 
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d N2 & = d N2 § 1 + e ( tr $Q ) ( tr $n ~ Q ) > ( 3 - 27 ) 



a=0 



£ 9k& k = (* " + 6^ n+fc ) . (3.28) 



fc>0 fc>0 



From the invariance of the partition function we obtain the so-called Schwinger- 
Dyson equations of the one-matrix model: 



/f> fc tr$ fc -A=0, (3.29) 



and 



^£(tr$ Q ) (tr <f> n ~ a )^ = ^-(^kg k tr $ n+fc ^ with n > 0. (3.30) 

These equations can be rewritten as the following constraint equations, acting on 
the partition function of the model: 



L n Z = for n > -1, (3.31) 



where: 



^0 = ^% fc ^- + iV 2 , (3.33) 

fc>0 ^ 

It is easy to check that these operators satisfy the Virasoro algebra: 
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[L n , L m ] = (n - m)L n+m for n,m> -1 



(3.35) 



and are therefore called Virasoro constraints. Since n,m > —1 the term pro- 
portional to the central charge does not appear. We have considered analytic 
transformations, since these are the only ones that can be expressed through the 
potential and its derivatives w.r.t. the coupling constants. 

Before we derive the loop equations we have to introduce the notion of the 
loop operator. The basic observables in quantum gravity are loops [14,17], since 
these are the only gauge-invariant quantities. We distinguish between two kinds of 
operators, microscopic and macroscopic loops. The lattice scaling operator has the 
form: (A i g/iV)tr$ n . The insertion of such an operator into the path integral creates 
a vertex with n spokes. It is easy to see that a hole with n edges is created in the 
dual surface, since the Euler characteristic has been changed to X = 2 — 2h — B, 
where B is the number of holes and B — 1 for one insertion. These operators are 
usually called 'microscopic' loops because they have only a few lattice spacings of 
length. They contain all the information about integrals over the surface of local 
operators. The generating functional of these loops is the so-called 'macroscopic' 
loop operator 



«o(0 = \^\ rfiv2$ex P (-£ trV w) tre ^ ( 3 - 36 ) 

The macroscopic loop operators correspond to extended boundaries on the surface 
and they can be expanded in terms of the lattice scaling operators. We will omit 
the expectation values in our notation: 



«o(0 = ^t re <* = £^«» = ££ M <">. ,3.37) 

n=0 ' n=0 

The moments = (Ag/iV) tr $ n can then be expressed through the free energy 
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of the model*: 



4 0, = -A||^ = A B and „<"> = -A||f . (3.38) 

We have used that the dependence of F on go is trivial F = —go/Ag + . . .. 

We observe, that from the Lq constraint, we obtain the discrete string equation. 
This means a relation between the parameter Ag and the operators tr $ fc : 



A s = ^E^A«(^ tr$fc )- ( 3 - 39 ) 

fc>0 ^ ' 

This equality can also be written as an integral equation, as will be done in section 
2.3. This will be important in order to classify the different multicritical behaviors. 

The L-i constraint takes the form: 

E^A B (^tr$ fc -A=0. (3.40) 
fc>i ^ ' 

This relation contains a non-trivial information if the potential of the model is 
arbitrary, and will be important in the next chapters. 

To find the solution of the Hermitian one-matrix model we can use the loop 
equation. It has been used by Kazakov [14] in the planar limit to analyse the 
critical regimes of (3.1); they provide us with non-perturbative information of the 
matrix model, as shown by David [63]. 

The loop equation can be determined from the Virasoro constraints. These 
constraints can be expressed as non-linear differential equations for the free energy. 
The L n constraint takes the form 



* Here we have rescaled F by a factor of TV 2 , i.e. we set Z — e N2p ; this is useful in order to 
guarantee a finite zeroth-order free energy. 
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Fig. 7: Geometrical representation of the loop equation. 



V" k gk dF + 3. Y ( d2F + at 2 — — ^ - (3 41) 

^ k dg k+n N 2 £-^\dg n - k dg k dg k dg n ^ k J 

In the limit N — > oo we obtain the equation: 




(3.42) 



This means that (3.42) is equivalent to: 

n 

E^4 n+fc) -E4 n ^4 fc) = o. (3.43) 

k>\ fc=0 
Adding up we obtain, from (3.43), after simple manipulations: 

I 

£ik — 1 p 

E ^^fWO = y dl'u Q (l - 0«o(0- (3-44) 
*>i o 

This equation was first derived by Kazakov [14] using more phenomenological ar- 
guments and is similar to the Migdal-Makeenko equations of multicolor QCD. The 
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loop equation has a very nice geometrical interpretation in terms of joining and 
splitting of loops, as shown in Fig. 7. It can be solved if we introduce the Laplace 
transformed loop operator: 



00 00 (it) 

«o(p) = I dle^uo(l) = J2 %T = V + • • • • (3 ' 45) 

J k=0 y P 



In terms of this operator the (non-planar) loop equation, that is equivalent to the 
Virasoro constraints, is: 



MP? - V'(p)u Q (p) + ^X(P,P) = Q(P)- (3-46) 

Here Q(p) is a polynomial in p whose explicit form is not necessary and we have 
introduced the two-loop operator: 



X(P, V) = X] k+i'k+i with Xfc,/ = o o , (3-47) 

and V'(p) = J2k>l kgkP k ~~ l - This equation has been solved by Kazakov in the 
naive iV — > 00 limit. For higher genus and even potentials the solution was found 
by Ambjorn and Makeenko [64], using the genus-expansion proposed by David [63]. 
We will now see how we can solve them on the sphere for an arbitrary potential 
[37] and how we can identify the multicritical points, the expressions for the scaling 
operators and the critical exponents. 
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3.3. Solution of the Planar Loop Equation for Arbitrary Poten- 
tials 

The solution of equation (3.46) in the limit N — > oo is 



The sign in (3.48) can be determined from the requirement that uo(p) is perturba- 
tive in the couplings and that its large p expansion contains only negative powers 
in p. If we choose one single cut dominating the scaling region, we can make the 
following ansatz 



In the above expression (x, y) are to auxiliary variables that will be expressed 
through the parameter A#. M(p) is a polynomial in p that will be determined in 
such a way that the loop equation is satisfied and that the loop operator has an 
expansion of order as follows from (3.45). 

The steps that we have to follow to determine the solution are: 

> determine x and y in terms of A#, 

> fix the expression for M(jp), 

> take the scaling limit, 

> compute the expressions for the scaling operators and the critical exponents, 
such as the string susceptibility. 

First we would like to derive an expression for the 'cosmological constant'* Ag, 
the so-called string equation. The loop operator uq(j>) has the following expansion 

* Actually, this is not the cosmological constant but it is related with it through (3.16); we 
will be sloppy with this terminology. 




(3.48) 



«o(p) = v '{p) - M(p)y/(p-x)(p-y). 



(3.49) 
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«o(p) = C(x, y) + + O (Jpj . (3.50) 

From (3.45) we obtain C(x,y) = 0. Dividing this expression by yj(p — x) (p — y) 
gives 



, "° = v " (p) -M(p) 



(3.51) 



+ 2A g (x,y) + (x + y)C(x,y) + Q / J_ 



which implies 



C(x, y) = -4> ^ - ^ = 0, (3.52) 



oo 

Equations (3.52) and (3.53) can be used, in principle, to rewrite x and y in terms 
of the single physical parameter Ag. These equations have appeared previously 
in the matrix model literature (see for instance [65] and references therein) in a 
slightly different form. 

Note that the piece proportional to (x + y)/2 in (3.53) vanishes when (3.52) 
is enforced. Then one recovers the expression for Ag, which is usually found in 
the literature. However, in order to compute partial derivatives it is essential to 
use the complete form (3.53). The above equations are equivalent to the planar 
discrete Virasoro constraints L_i and Lq respectively. 
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We now proceed to the computation of M(p); it can be expanded in powers of 

P 



M(p) = J2M nP n . (3.54) 



n>0 



From (3.51) we obtain explicit expressions for the expansion coefficients: 



M„ = -<^ ,'- V W , (3.55) 



so that we obtain the integral representation 



M(p) = -<f£- ; w (3-56) 



2« ( g _ p)v/(g-a;)(g-|/) 



Therefore, we have expressed uq as a function of Ag only, since M(p) is a function 
of (x, y) through (3.56) and (x,y) are functions of Ag, since (3.52) and (3.53) hold. 

In order to derive the basic observables of our theory, the scaling operators, 
we would like to obtain an expression for the derivative of uq(j>) with respect to 
the cosmological constant Ag. Under a variation (dx, dy) compatible with the 
constraint (3.52) duQ is given by 



du Q (p) = -L (-±M(p)dA - AdM(p) S j = d ~Y + o[^). (3.57) 



We have introduced the notation A = (p — x){p — y). Since the p-dependence of 
the bracket cancels out, we obtain the expression 
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dA B = -]-M(p)dA - AdM(p). (3.58) 



Thus, the result for uq is: 

du (p) 



dA B y/(p-x)(p-y) 
We are now in the position to analyse the theory in the scaling limit. 



(3.59) 



First we would like to see how this limit could be taken if the potential is 
arbitrary. For this purpose we use an alternative definition of the macroscopic 
loop operator (3.37): 

« (0 = lim -^4 n) > (3-60) 

where / = na?/ m is the renormalized length. Here m is a positive integer related 
to the order of criticality as we will see and a is the cutoff. As n goes to infinity 
the lattice spacing a will go to zero, but I has to be kept fixed. This operator has 
the same continuum limit as the lattice loop operator (3.36) up to a multiplicative 
renormalization. We will therefore use the same notation for both of them. The 
continuum limit has to be defined in such a way that (3.60) makes sense, since this 
operator is the basic geometric observable in 2D-quantum gravity. This will be our 
guiding principle. Equations (3.59) and (3.60) imply 



d AB u (l) = hm -L^ / dpP " (3.61) 

n^oo 2m A B J y/(p-x)(p-y) 

c 

where the contour C encloses the cut (x,y). We may assume \y\ < \x\ < 1 without 
loss of generality. Taking C as the unit circle around the origin with p = e t9 we 
have the identity 



•57 



1 / dp p n If ^ e%9 ex P 2//m ^) 



2lii J ^{P~ x){p-y) 2ttJ ^1 ( e i0 — x ){e ie — y) 
c c 



Changing to the new variable z = a 2 / m 6, eq. (3.62) becomes 



(3.62) 



dze ilz 

<W0 - ■'- ' 



(exp (ia 2 / m z) — x)(exp (ia 2 l m z) — y) 

(3.63) 

Here we have taken x c — 1 without loss of generality. This expression will vanish 
unless the integral itself is of order a" 1 /™, i.e. if x approaches the critical value 

x c = 1 as a 2 / m 



x = 1 - a 2 l m u + . (3.64) 

The variable u + is usually called the 'specific heat'. We will later see how it 
is related to the free energy of the model. Since x and y are not independent 
variables, y will approach its critical value y c (\y c \ < 1) at the same time 



y = y c + a 2 / n u-. (3.65) 

Here n is another positive integer that may be different from m. The concrete value 
of y c is fixed by x c through eq. (3.52). The previous two equations are important, 
since they indicate to us how we have to take the continuum limit. 

The integral in (3.63) is dominated by the region p ~ x c — 1, where the contour 
can be deformed into a straight line*: 

★ Here we have absorbed a factor y/1 — y c into the definition of k. 
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+ioo 

1 / e lz 

ndtunil) = / dz — — . (3.66) 

uv ; 2m J y/z + u+ v ; 



-too 



Here we recognize the definition of the inverse Laplace transform 

u (z) = kC[u {1)], (3.67) 

where 

d t u (z) = - 1 . (3.68) 

V z + u + 

The previous result can of course also be obtained from (3.59) if we take into 
account the scaling law 



p = 1 + a 2/m z (3.69) 

together with (3.64) and (3.65). The loop operator in the scaling limit uq(z) is 
connected with the discrete one through the relation: 



u (z) =a~ 2+1 / m u (p). (3.70) 

There are therefore two different ways of viewing uq(z): as the Laplace transformed 
of the macroscopic loop uq(1) or as the continuum limit of the loop operator uq{j>). 

From (3.66), we get the expression for the macroscopic loop after integration 



1 p -lu+ 

d t u (l) = TT . (3.71) 

Note that dfUo(l) is independent of the variable U-, so that u + is the 'physical' 
scaling variable . This result holds as long as \y c \ < \x c \ independently of the values 
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of m and n. It is easy to see that for \y c \ = \x c \, the dominant end-point will be 
the one with the highest order of criticality. 

We would like to obtain the expressions for the scaling operators and the form 
of the string equation. In analogy to the discrete partition function (3.1) we can 
describe the non-perturbative partition function of a general massive model in this 
limit in a formal way as follows: 



^(*0 >*l ' • • •) — 




(3.72) 



Here a* are the scaling operators in the continuum theory, while tj' are the renor- 
malized coupling constants. Therefore we can derive the correlation functions of 
the scaling operators in the background t% = b n , for n > 0, from logZ as 
follows 

K ■ ■■<) = ^ T^T ^Z(t+,t+, . . .)\ lt=K . (3.73) 

Ul ni . . . Ol ns 

The renormalized cosmo logical constant t has to be identified with tj". 

We first compute the one-point functions. The macroscopic loop operator 




(3.74) 



has an expansion in scaling operators (a£) [13,14,17] dual to the renormalized 
couplings tf that can be used to compute explicit expressions for these operators: 



00 ik i 

5o(0 = «D-)* +1 ifci< 3? ?>-^ ( 3 - 75 ) 

This formula is the analog of the formula (3.37) in the continuum. For the scaling 
operators we obtain the expressions 
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From the above formula it becomes clear that the specific heat is the two-point 
function of the puncture operator 3q~: 



= <W> = ^f- (3-77) 



Next we would like to obtain the form of the string equation for generic po- 
tentials. To do so, we have to obtain the corresponding equations to (3.52) and 
(3.53) in the scaling limit. Comparing the first derivatives of A B (x, y) and C(x,y) 
we find that they are not independent. Instead, 



d x A B = ^(x- y)d x C 



d y A B = ^(y-x)d y C. 



Similarly, we have the following identities 



(x - y)dl y k B = ~\{d x A B - d y A B ) 



2 



(3.78) 



(3.79) 



( x - y )dl y C = -{d x C-d y C). 
Scaling the first equation of (3.79) according to (3.64) and (3.65) gives 

d\ _t ~ (a 2 ' n d + t + a 2 ' m d^t), (3.80) 
with an identical expression for C. For a — > the r.h.s. vanishes and we find 
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dl_t = and d\_c = 0. (3.81) 

Here we have defined d± = d/du± and c = — a~ 2 (l — y c )C/2. Moreover (3.78) 
implies 



<9 + c = <9+i and d-c = -d-t. (3.82) 

One solution of (3.82) is given, if the renormalized cosmological constant t and c 
are expanded in terms of scaling operators 

t = -^J2?p u + + ^p u -) > 

(3.83) 

Equations (3.83) are the continuum version of (3.52) and (3.53) or equivalently 
of (3.39) and (3.40). A (more rigorous) prove of these equations can be obtained 
using the continuum loop equations. This will be done in section 5.2. for the super- 
symmetric case. The purely bosonic model can be obtained setting the fermionic 
coupling constants to zero. 

Adding and subtracting the equations in (3.83) we obtain the string equation 
in terms of the physical variable u+, which is the one that appears in the loop 
operator (3.71) 



t = -EW- ( 3 - 84 ) 

Therefore, we have a parametrization of the string equation and the loop operator 
in terms of the '+' variables. The reason for this is that for \y\ < \x\ the continuum 
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limit is controlled by m. The converse is of course true for \y\ > \x\ and we would 
have a dependence of the loop operator and the string equation on the '— ' variables. 
From now on we will omit the subindices in our formulas. 

At the m th multicritical point the planar string equation has the following form 



u m 



t n u n . (3.85) 



n>m 



The cosmological constant therefore satisfies 



dt 
du 



u=0 



d 
du 



TO— 1 



t 

m—1 



= 0. 



(3.86) 



u=0 



Using the above string equation it is easy to see, that the string susceptibility 
has the value T s t r (0) = —1/m. The scaling dimension of an operator (a^) is given 
by dfc = k/m. These values can be compared with those obtained from the Liouville 
approach in chapter 2. They correspond to the (2m — 1, 2) minimal CFTs coupled 
to 2D gravity. 



3.4. Basic Concepts on Integrable Hierarchies 



Two-dimensional quantum gravity can be formulated non-perturbatively if we in- 
troduce the double scaling limit [13], as mentioned in section 2.1. It turns out that 
in this limit a close connection to integrable hierarchies of the KdV type appears 
since, for example, correlation functions of microscopic operators can be expressed 
in terms of KdV flows. 

In this section we are going to explain some basic notions on integrable hier- 
archies and the relation to two-dimensional quantum gravity will be the subject of 
the next section. A more complete description of the field is given in some good 
reviews on integrable hierarchies, as for example [66,67,68,55]. 
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The following definition of an integrable system dates back to Liouville: A 
Hamiltonian system with a 2iV-dimensional phase space is integrable if and only 
if there exist exactly N functionally independent conserved quantities that are in 
involution (that is, the Poisson brackets of these conserved quantities with one 
another vanish). If the system under consideration is infinite-dimensional this 
would imply the existence of an infinite number of conserved quantities. 

The KdV equation is one example of such a system. It is a non-linear equation 
in one space and one time coordinate (x,ti): 



du du (3 87) 

dt\ dx dx 3 ' 

Here u = u(x, t\) is the dynamical variable and k is a free parameter that we have 
introduced for later convenience. This equation can be written as an Hamiltonian 
equation: 



Ou 

— = {«,Wi} 2 = {«,W 2 }i. (3.88) 

Oil 

Here Hi and H2 are the Hamiltonians of the system. It is a peculiarity of the KdV 
equation that we can find two Hamiltonians and two Poisson brackets, so that 
(3.88) is equivalent to (3.87). The two possible definitions of the above Poisson 
brackets are: 



{u(x),u(y)}i = d x 5(x - y) 

(3.89) 

{u(x),u(y)} 2 = (k 2 x + 2u8 x + 2d x u) 5(x - y). 

They satisfy the Jacobi identities and are antisymmetric. The second structure is 
identical to the classical version of the Virasoro algebra. The corresponding two 
Hamiltonians are 
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Hi = J dx , 

— oo 

(3.90) 

— oo ^ ' 

respectively. 

The solution of the KdV equation is unique for a given initial condition. It can 
be shown that waves without dispersion, i.e. which maintain their shape as they 
move, are a class of solutions. They are called solitons, and imply the existence of 
an infinite number of conserved quantities TC n . These conserved quantities satisfy: 

^ = {H n , H 2 }i = {H n , Hih = 0. (3.91) 
ati 

The simplest example is of such a constant of motion is 



00 




(3.92) 



— 00 



Using the KdV equation (3.87), it is easy to prove that other constants of motion 
are for example Hi and H.2- One can explicitly check that these three quantities 
satisfy the recursion relation 



U 2 8l + 2ud x + 2d x u) 5 -^- = dj-^ for n = 0, 1, 2, (3.93) 
v ' bu bu 



where H-i = and the functional derivative is defined as 
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5F faffl = l [ ™ ~e (F[M(;r)] + * ~ y)] ~ F ^ x ^ ■ ( 3 - 94 ) 

The first three conserved quantities can be obtained by taking into account 
simple symmetry principles. We can identify 7i\ with the momentum since it 
generates space translations with respect to the first Poisson bracket 



du 

{u > Hl}l = fa- (3 ' 95) 

Similarly H2 is the generator of time translations and can therefore be identified 
with the energy 



{^2)1^. (3.96) 

It can be proven inductively that we can choose all conserved quantities 7i n 
(3.91), so that (3.93) holds for n£N. Furthermore this recursion relation guar- 
antees that these quantities are in involution: 



{H n , H m }! = {H n , H m } 2 = 0. (3.97) 

The relation (3.93) is usually written as a recursion relation for the Gel'fand-Dikii 
polynomials Rk[u] that are closely related to the constants of motion: 



{K 2 dl + 2ud x + 2d x u))R k ^[u} = d x R k [u] with R k [ u } = 6 -^±, (3.98) 



with Rq — 1/2, Ri — u, R2 — Su 2 + n 2 d 2 u, .... Starting with the expression for 
Ho in eq. (3.92) we can use (3.98) to determine all the other charges. 
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Each of the conserved quantities can be considered as an Hamiltonian that 
generates an evolution equation, that is a generalization of (3.87) : 



du 
dt k 



{u,H k }2 = {u,Ti k+ i}i, 



(3.99) 



or equivalently 



du 
dt k 



d x R k+1 [u\. 



(3.100) 



The function u can be thought of as depending on an infinite number of time 
variables u = u(x, t\, . . .), where we define x = to- These are known as higher flows 
of the KdV hierarchy. 

Before we establish the connection between this formalism and two-dimensional 
gravity we would like to introduce a suitable formulation that can be used for 
all the (p, g)-minimal models coupled to two-dimensional gravity: the Lax-pair 
formulation. The formalism of pseudo-differential operators that we now explain 
will be relevant to describe the supersymmetric extensions of the KP hierarchy in 
section 4.6. 

The p th -reduction of the generalized KdV hierarchy is defined in terms of a 
differential operator: 



Here U{ are a priori arbitrary functions depending on the time variables (to, ti, t2, ■ ■ ■) 
We can consider the flows 



★ It is always possible to set the order (p — 1) to zero by conjugation. 



p-2 




(3.101) 



i=0 
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^ = [H n ,L], (3.102) 

generated by the infinite set of commuting Hamiltonians H n . The Hamiltonians 
are represented as fractional powers of the positive part of L n 



H n = (L n /P)+- (3.103) 

A fractional power of L is defined as a Laurent series in the differential operator 
d x . These are in general n-th pseudo-differential operators, i.e. operators of the 
form 



0= J2 ( 3 - 104 ) 

i=— oo 

where the operator d~ l acts on a function / as follows: 



i=0 

The part of the sum that contains only positive powers in d x is denoted by 0+: 



0+ = J2°i d x ( 3 - 106 ) 



i=0 

and similarly O- = O — O+. The residue of a pseudo-differential operator is 
defined as the coefficient of d~ x in the Laurent expansion of O 

res £> = £>_!. (3.107) 
If n is a multiple of p, the flow is trivial 
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— = for n = pq, 
dt n 



(3.108) 



since H n = L q commutes with L. 

The simplest reduction is the 2-reduction where 

L = dl + u, (3.109) 

whose first non-trivial flow is the KdV equation (3.87). This operator plays an 
important role in the non-perturbative description of the one-matrix model as we 
will see later. The flows computed at the beginning of this section can now be 
obtained from (3.102) comparing the powers of d x of the right- and left hand side 
of this equation. 

The generalized KdV hierarchies are special cases of the KP hierarchy, which 
is defined in terms of an arbitrary pseudo-differential operator 

oo 

Q = d x + Y,Q^ 1 - ( 3 - 110 ) 

i=i 

The KP hierarchy is defined as the infinite set of commuting flows on the phase 
space of the operator Q 

^ = m,Q] = -m,Q\- (3.1H) 

The p th KdV hierarchy satisfies the condition that Q p = L is a pure differential 
operator, i.e. Q p _ = 0. 

The absence of a term Q$d® in Q is sufficient for the existence of a pseudo- 
differential operator 
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S = 1 + ^ Sid' 1 with Q = SdxS' 1 . (3.112) 
i=i 

In terms of this operator the flow equations become 



— = -[SdZS-^S. (3.113) 

To the p-reduction of the generalized KdV hierarchy (3.102) we can associate 
the so-called r-function that is related to the differential operator L as follows: 



resi " /p = aS; logT - (3 - 114) 

The r-function completely determines the operator L since for % — 1, . . . ,p — 1 we 
have 



res L*/p = % -u v -i- X . (3.115) 

If we consider the KdV operator, the variable u, as a function of which the 
flows are formulated, is completely determined by the r-function: 

u = 2d 2 x log r. (3.116) 
The r-function of the KdV hierarchy satisfies the recurrence relation 



d 2 x - logr= ( K 2 d 3 x + 2ud x + 2d x u)d x —\ogr. (3.117) 

Otn+l Ot n 
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3.5. Double Scaling Limit, KdV and Continuum Virasoro Con- 
straints 

The r-function is useful to describe the partition function in the matrix-model 
approach to 2D quantum gravity as we will now see. The important result that 
provides the link between the non-perturbative approach to 2D-quantum gravity 
[13,17] and integrable systems is that the partition function of the one-matrix 
model is the square root of the tau-function of the KdV hierarchy [16]: 



Z(tQ,t!,...) =r 2 (t , ti,...). (3.118) 

Since u is related to the r-function as stated in (3.116) it can be identified with 
the 'specific heat' or equivalently the two-point function of the puncture operator 
(3.77). Recall that this means that the specific heat satisfies the KdV flow equations 



^ = DR n+1 [u}. (3.119) 



Here k is the renormalized string coupling constant and D = d/dto. 

The insertion of a n corresponding to the differentiation with respect to the 
coupling t n is then identified with the n th KdV flow of u: 



(a n a a ) = — = DR n+1 [u}. (3.120) 
dt n 

Therefore, the recursion relations of the i2fc[u]'s imply recursion relations for scaling 
operators. The non-perturbative string equation can be written in terms of the 
Gel'fand-Dikii polynomials 
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t = -J2(2k + l)t k R k [u]. 
k>l 



(3.121) 



This equation is a generalization of (3.84) to higher genus. At the m -multicritical 
point all t k vanish except t m . This coupling constant can be normalized in such a 
way that for the first values in m the string equation takes the form 



m — 1 
m = 2 



t = u 

K 2 

t = u 2 H d 2 u. 



(3.122) 



The second equation is known as the Painleve-I equation; it is the string equation 
of pure gravity. For k — (planar limit) we know that the specific heat behaves in 
the pure gravity case as u(t) = —\fi\ this is thus the initial condition for the genus 
expansion of the solution. The first terms of this expansion are 



u{t) = -Vi+ -t- 2 + —K 2 t- g / 2 + .... (3.123) 
w 24 1152 v ; 

It turns out that the string equation can be written as the following constraint 
equation on the square root of the partition function, i.e. on the r-function [18,19] 




(3.124) 



Using the relation between the partition function and the r-function, one can show 
that the partition function satisfies the constraints [18,19] 
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L n r = for n > -1, (3.125) 
with the L n operators defined as 

m=l 



L n = ]T m +- t m - + -^^ s . (3.128) 

m=0 ^ _/ m+n ™-i 



t'm+n 2 dt m —idt n 
These operators satisfy the Virasoro algebra: 



[L n , L m ] = (n - m)L n+m for n,m>-l, (3.129) 

so that they are the Virasoro constraints in the double scaling limit. They are the 
continuum version of the constraints (3.31). Notice that for the discrete theory 
the constraints acted on the partition function itself and not on the square root. 
The conditions (3.125) are due to the fact that the the recursion relation for the 
T-function (3.117) implies: 




(3.130) 



Together with the string equation (3.124) these equations imply the constraints 
(3.125). 



73 



The expressions (3.126), (3.127) and (3.128) can be written as the Virasoro 
generators of a ^-twisted free scalar field x(z), i.e. with x(e 27Tl z) = —x(z). The 
mode-expansion of x(z) is 

dx{z)=Y. a n+\z- n - 3 >- (3-131) 
Since the energy-momentum tensor is given by 



T(z) =: \dx{zf : + J-^ = £ L n z~ n ~ 2 , (3.132) 

we obtain the following correspondences between the modes of the free boson and 
the coupling constants of the model: 



— [n-\ — J — and a n +- = K ~0~ ^ or n>0. (3.133) 
2, J K 2 otn 



The Virasoro constraints in the double scaling limit can be shown to be equiv- 
alent to the continuum non-perturbative loop equation of the one-matrix model 
[18,19] in a way similar to what was already done in the discrete theory in section 
3.2. This equation has the form 



Uq(z) + k?x(z, z) = polynomial (z), (3.134) 
where we have introduced the two-loop operator 



1 

X(z, z) = ^ + 2^ 2 J2 52(Wn-k) (3-135) 

n>0 fc=0 
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and the loop operator takes the form 

«o(*) = ( n + I) + k2 H *~ n_l K>> (3.136) 

n>0 ^ ' n>0 

where uq(z) and x\ z i z ) are the continuum limit of the discrete loop and two-loop 
operators respectively: 



«o(p) = «o(p) - 



Hp) 



2-- 



: w (^) and j^x{p,p) 



2 4-^ 
/c a m 



x(^)- 

(3.137) 
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4. THE SUPERSYMMETRIC EIGENVALUE 
MODEL FOR GENUS ZERO 

In this chapter we begin to analyse how the previous formalism could be generalized 
to the simplest N — 1 SCFTs coupled to 2D-supergravity. For this generalization 
the Virasoro constraints satisfied by the partition function will play a central role. 
We start with the construction of the supersymmetric eigenvalue model found by 
Alvarez-Gaume et al. [30]. 

4.1. Construction of the N — 1 Supersymmetric Eigenvalue Model 

We will take as a guiding principle the one-matrix model formulated in terms of 
eigenvalues and construct, in analogy, the N = 1 supersymmetric generalization. 

4.1.1. The One-Matrix Model in Terms of Eigenvalues as Guiding Principle! 

In analogy to the continuum constraints, we can write the discrete Virasoro gener- 
ators (3.31) as the modes of the energy-momentum tensor of a free boson dx. We 
can introduce an infinite set of creation and annihilation operators: 

that are the modes of dx{p) 
tensor, 

T{p) = l - : dx(p)dx(p) := f^p^- 2 L n , (4.2) 
are then represented by the operators (3.31). The discrete-loop equation (3.46) 



N 

.ng n tor n > (J 



A B V2 

(4.1) 

A B V2 d 



for n > 0, 



N dg n 

^ nGz CinV n l - The modes of the energy-momentum| 



7G 



becomes: 



Z 1 T(p)Z = polynomial (p). 



(4.3) 



On the other hand, we can start with the free boson with anti-periodic bound- 
ary conditions and show that we can determine the partition function of the one 
matrix model. First, we have to construct the potential of the model in terms 
of eigenvalues and then determine the measure. It is easy to see from (4.1) that 
precisely the potential is related to the negative modes of the free boson 



The loop operator can be determined from the positive part of the mede ex- 
pansion of x 



Next, we can write the partition function (3.1) in terms of the eigenvalues of the 
$-matrix (Ai, A2, . . . , Ajv), but leaving the measure A 2 (A) undetermined: 




(4.4) 



fc>0 




(4.5) 




(4.6) 



The constraints (3.31) yield a differential equation satisfied by A 




(4.7) 
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whose solution up to a constant is the expected Vandermonde determinant 

A = II( A *- A i)- ( 4 - 8 ) 

i<j 

Therefore, if the potential of the model in terms of eigenvalues is known, we can 
use the Virasoro constraints to determine the measure. 

4.1.2. The Generalization to the N = 1 Supersymmetric Eigenvalue 
Model 

In the supersymmetric case we proceed by analogy with the above arguments 
to obtain the expression for the partition function in terms of eigenvalues and 
the corresponding loop equations. We begin with the derivation of the discrete 
superloop equations. This will be done by taking a c = 1 free massless superfield 
X(p, n) = x(p) + Hip(p) as a generalization of the free boson. These fields have a 
mode expansion 

dx(p) = ^2 a nP n l and ip{p) = ^ b r p^ r ~^, (4.9) 

where the fermion is in the Neveu-Schwarz sector. The energy-momentum tensor 
is 

T(p, n) oc DXd p X = i)d p x + IT : (d p xd p x + dpipip) :, (4.10) 

where D = d/dU + Ud p is the superderivative and (p, II) are the two coordinates of 
the superspace. The loop operator w(l, 6) in the supersymmetric theory depends on 
two variables (I, 8) that are the even and odd length respectively. They characterize 
the boundary of a superdisk. We can define the Laplace transformed loop operator: 
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it; 



oo 

(p, n) = v(p) + Uu(p) = J dl J d9e- pl - m w(l, 0). (4.1i; 



The operator w(p,U) can be identified with the field X+(p, II), so that it has an 
expansion in negative powers of p: 



7} 7) 



The analogy with the bosonic theory suggests that we can identify the bosonic and 
the fermionic modes as bosonic and fermionic couplings respectively 



A $ 9 N f mo 

ap = ~Jl~dip ®-p = -j^P9p, for P = 0,l,2, 



(4.13) 

If we write F = A^~ 2 log Z then the moments and can be identified with 
derivatives of the free energy 



u^ = A s „(») = -A||^ V (») = _A|^1. (4.14) 

C<7n W§ 



As a basic postulate we take the (non-planar) discrete superloop equations to be 
described by the above energy-momentum tensor [30] 
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Z _1 T(p, U)Z = polynomial (p). 



(4.15) 



In terms of u(p) and v (p) this is equivalent to a system of two equations: 

Hp) - v\p)f + (v( P ) - t(pMv(p) - £(p)) + ^ (x BI W) + x ff (p,p)) = Qo, 

(4.16) 



u(p)v(p) - V'(p)v(p) - Z{p)u(p) + ^ x BF (p,p) = Q x . (4.17) 
Here we have introduced the potentials: 



V{p) = 9kP k and t(p) = &+1/2P* (4-18) 

fc>0 fc>0 

and the two-loop operators are defined as: 



X M= ^ jfi+lgl+l dgkdgi -- 



k,l>0 ■ 



BF/ \ d 2 F 



x ff (p.p) = EE^i( 



A.% rn \ d 2 F 

r 



n>l r= 



The quantities Qo and Q\ are polynomials in p. Although their explicit form can 
be computed they will not be needed. In terms of the original loop operator w(l,9) 
eqs. (4.16) and (4.17) take a form similar to (3.44): 
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VlCw(l, 9) + 2K,Vw(l, 9) = (wo Vw)(l, 9), (4.20) 



with 



The convolution between two superfunctions fi(z, 9) and f2(z, 9) is defined accord- 
ing to: 

2 

(fog)(z : 9) = J d9' J f(z',9')g(z-z',9-9')dz'. (4.22) 



Now we would like to see what the discrete super- Virasoro constraints and the 
partition function on which they act look like. The analogue of the potential of the 
one-matrix model is identified with the negative mode expansion of the superfield 
DV(p,U) ~ DX- 

v( P , n) = J2(9kP k + 4+in/). (4.23) 

fc>0 

The previous formula suggests the introduction of a superpotential that in terms 
of the even eigenvalues Aj and the odd eigenvalues 9{ takes the form 



N 

n\ «) = E Ek* A ? + tk+^h (4.24) 

fc>0 z=l 

Writing the partition function as: 
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Z = / ft d\id9iA(X, 9) exp [~^ v (^ > ( 4 - 25 ) 

we can determine the explicit form of the measure A(A, 6) by imposing the super- 
Virasoro constraints. The explicit representation of the super- Virasoro operators 
using (4.13) as differential operators is 



d A, 9 A 2 ™ «9 2 



fro 2<9 ^ fro ^ fro 



2N 2 ^ V2 / d&dtn-r 



for n > — 1 



(4.26) 



These operators satisfy a closed subalgebra of the N = 1 superconformal algebra 
in the Neveu-Schwarz sector. The partition function of the supersymmetric model 
is annihilated by L n for n > —1 and by G n _i for n > 0. Since the last equation 
in (1.48) holds, it suffices to implement the G n _i constraint. This leads to a set 
of equations of the form: 

whose unique solution, up to a multiplicative constant, is: 
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A(\,8) = Y[(Xi-\ j -6 i e j ). (4.28) 

i<j 



Hence the model we would like to solve is: 



/N / oo oo \ 

J] dXidOi - Xj - BiBj) exp - — E + ^+|^ A • 

i=l i<? V 5 fc=0 i=0 / 

(4.29) 

Finally, we would like to remark that the loop operator can be explicitly written 
in terms of the eigenvalues as: 



™M)^I> /A<+ ^ (4-30) 

i 

In the next section we begin the analysis of this model in the large- N limit. 

4.2. Solution to the Planar Superloop Equations: General Poten- 
tial 

We are going to calculate the expressions for the bosonic and fermionic loop op- 
erators that follow from the previous model in the spherical topology. The planar 
loop equations follow from (4.16) and (4.17) for iV — > oo: 



(u(p) - V'{p)) 2 + (v(p) - t(p))' (v(p) - tip)) = Qoip), (4.31) 



Hp) - e(p)) Hp) - v'(p)) = gx(p). (4.32) 

As in the bosonic case we look for the one-cut solution. We will consider a general 
bosonic part of the superpotential so that the one-cut ansatz takes the form: 
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u(p) = u (p) + u 2 {p) = V'ip) - M(p)VA - -gl, (4.33) 



v(p) = £(p) - (4.34) 

with A = (p — x){p — y). Here uq{j>) coincides the solution of the purely bosonic 
model (3.59) and determines the planar solution in order zero of the fermionic 
couplings. The subindex in (4.33) indicates the order in fermionic couplings. As 
we will now see, the superloop equations can be solved taking for the loop operator 
an expansion that is at most quadratic in fermions. It has been argued in [30] that 
the solution to the superloop equations is unique, so that this property is inherent 
in the model. It will also appear non-perturbatively [40], as we will see later, 
and means that correlation functions involving more than two fermionic scaling 
operators will vanish. 

To obtain the complete solution of the model we have to determine the func- 
tions M(p), A(p) and N(p) in (4.33) and (4.34); M(p) has already been determined 
in (3.56). Here it will be written in a simpler form, in terms of recurrence relations. 
We introduce the notation 



x = S+VR and y = S-VR. (4.35) 

To write down the explicit form of M, N and A, we note that any analytic function 
f(p) can be written in the form: 



f(p)=f-(A) + A'f + (A) with A ' = ^ = 2(p-S), (4.36) 

where we have split the function f(p) into two terms of opposite parity with respect 
to the change (p — S) — > (S — p). More concretely 
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M(p) = M _ (A) + A'M+(A), 

N(p) = AN~ (A) + A'iV + (A), (4.37) 
A(p) =A-(A)+A'A + (A). 



We can expand M, iV and A in powers of A 



M±(A)=J>±A fc , A/± = J>±A fe , A±(A) = J>±A*, (4.3? 

fc>0 fc>0 fc>0 



and we have to calculate the expansion coefficients. To determine A we substitute 
the above expressions in (4.33) and (4.34) and require that the left-hand side of 
(4.31) and (4.32) be polynomials in p. After some computations we obtain that 
A~(A) and A + (A) are completely determined by and Oq . The results are 



2R 

A- = a = -—— 2 —, + , 2 K n o"o - 4:Rm+n+n+), (4.39) 

(m ) £ — AR(rriQ ) z 



2R 

A+ = a+ = -—— 2 —, + (m Q n+n+ - m+n Q n+). (4.40) 

(m ) A — 4R(m,Q Y 



Therefore, once we have determined the coefficients and we will have the 
complete solution to the planar model. We split M(p) from (3.56) in M ± (A), and 
find: 
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M-(A) = & — ^ t{1 - tS) (4 41) 

f 2m ^(i _ tx)(l - ty) (1 - te)(l - ty) - t*A { ^ 



o 



M+f A) = I I dt 9{t ~ 1] - t (4 42) 

1 ; 2 7 2th ^(i - - ty) (1 - tx)(l - ty) - t 2 A- v ' ; 







Here we have introduced g(p) = pV'(p) and transformed the integrals so that the 

4 



contour goes around zero. Expanding in powers of A we obtain for m¥ explicit 



formulae that satisfy the following recursion relations: 



dm k {, 3\ ± dm+ ( 3\ _ (A ao . 



dm~ ( 3 \ 

= 4(fc + l)m+ + AR\k + - J m+ +l . (4.44) 



Therefore, we only have to determine the initial values m Q , and the values of 
the rest of the m^'s follow from the previous relations. This can be done by 
differentiating (3.53) with respect to x and y and comparing to (3.56). This gives: 



d x K s = \{x-y)M{x) and d y A s = \{y - x)M{y). (4.45) 
On the other hand by (4.38): 

M(x) = rriQ + (x — y)rriQ, and M(y) = ttIq + (y — x)m,Q. (4.46) 

This gives in terms of the single function Ag(x,y) 
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2 2 
m o = jjrz -yj( d x~ d y) A s and m o = ^ x _ y) 2 ^ + d y) A s- (4.47) 

Therefore, M{p) has been determined as a function of As- 

We continue now with the calculation of N±(A) that goes along the same lines. 
We have the following representations as contour integrals: 



r dt tit- 1 ) l-st 



v ; 2J2ni J (1-tx )(l-ty )-t 2 . 



(4.49) 



o 



(1 -ta)(l -ty) -t 2 A' 
Expanding in powers of A we obtain that the n^'s satisfy the recursion relations: 

' k + 7;)n k+v (k + - )n k , 



dR - V ^ 2 r k+v ds - V 2 ^ 



OR 



(* + <+i> = ^ k + + 4i? (* + I) 4+2- 

(4.50) 



Here we see that all the n^'s can be written in terms of Uq (the fermionic scaling 
variable), which can be defined through a contour integral in the same way as the 
cosmological constant in (3.53): 



+ _ i / dt at- 1 ) 





,+ 



n+ = - <f> , (4.51' 

o 2/ 2mty/(t-x)(t-y) 



The function has a property that will later be relevant in the continuum limit 
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(a: - y)d xy n+ = ^(d x n£ - d y n%). (4.52) 

The complete solution of the planar model has now been computed. 

We can use these results to obtain explicit expressions for the loop operators 
in the discrete theory. The above recursion relations and some amount of algebra 
determine the total derivative of the fermionic loop operator v(p) with respect to 
R: 



dv(p) I (dn+ dSdn+ 2RdS \ n+A> 

The order zero in fermionic couplings of the bosonic loop operator uq(p) has already 
been determined and is given by (3.59). The expression for U2{p) follows from 
considerations similar to those for the fermionic loop operator: 



u 2 {p) 



4R 



(m ) 2 - AR(m^) 2 



m, 



. dn,Q 



4Rm+ dn ° 



'"0 



dR 



A' 



in. 



driQ 
~dR 



m, 



+ dn o 
OS 

(4.54) 



n 



+ 



— ^— I 

A 3 /2 



In order to take the continuum limit we want to transform the derivatives with 
respect to R and S appearing in U2(p) into derivatives with respect to A. Using 
implicit derivatives we obtain 



dx 



R- 1 ' 2 



OA m Q + 2^m^ 



dy_ 
OA 



-2y/R, 



m 



+ 



(4.55) 



and since 
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d i f d d\ d i(d d . . 

and ^5 = —T. I in ~ IT.) ( 4 - 56 ) 



dR x — y \dx dy J dR x — y \dx dy 



we obtain the result 



Since the continuum limit of uq(p) has been computed in (3.68), we now determine 
the continuum limits of U2(p) and v{p). All we have to do for that is to determine 
the behavior of n+ in the scaling limit. We denote by r the scaling function 
corresponding to Uq . To determine the scaling of we notice that the scaling of 
U2(p) should be 



u 2 (z) = a 2+ ™u 2 (p) (4.58) 
as determined from the scaling of uo(p) in (3.70). Therefore Uq should behave as: 



4 = (g^) (4.59) 

where the sign is fixed by the continuum limit of the fermionic-loop operator. We 
notice that eq. (4.52) gives, in the continuum limit, an equation for r: 



d 2 r 

,d+-r = 0. (4.60) 



du + du- 

A solution to this equation is given by 
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r = ro + E( ?r M + ^ fc )- (4-61) 
fc>0 

Here = are the renormalized fermionic couplings. For simplicity can we set 
consistently r k = 0. We will omit the + and — signs introduced for the specific 
heat in section 3.3. The complete bosonic-loop operator is then: 



du(z) 1 d ( rdtr , 

In a similar way we get the continuum limit of the fermionic loop operator. Using 
(4.53) and transforming the derivatives, we get in the scaling limit: 



^ = -r^=^A where v(z) = a-^JT^vip). 

at I y/z + u at J 

(4.63) 

The factor a 1//m can be absorbed into the bare fermionic superlength 6 B 



w(l,6 B ) = u(l) + 6 B v(l), (4.64) 
because we can define renormalized superlength 



= 9 B a 1/m , (4.65) 

that behaves dimensionally as (length) 1 / 2 . 

Now we would like to obtain the expressions for the macroscopic loop operators 
from which we can obtain the bosonic and fermionic scaling operators. They are 
defined in a similar way as in (3.60): 
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u(i) = Mi) + Mi) = ^(E A * ) = J™o x^ (n) ' ( 4 - 66 ) 



u(Z) = lim (irOiX?) = lim T^ (n) - ( 4 -67) 

w ra->oo\^— j' */ moo Ag 

We can calculate the asymptotic behaviors using some standard formulas from the 
theory of Legendre and Gegenbauer polynomials [69,70]. To give an example, we 
start with the derivation of the uq operator, which has already been done in section 
3.3 through a contour integral. We can expand in terms of Legendre polynomials 



jMg) - v J_^l - Vf™W 2 P ( Z±A\ J_ (a eg) 

OAs ~ OAs ~ ^ XV) Pn { } 



n>0^ ~ AS n>0 

so that the expansion coefficients become 



= W n/2 Pn ( Z±±) ■ (4.69) 
OAs 1 y> \2^xyJ 1 ; 

Using the relation to the hypergeometric function F 

P n (z) = F(-n, n + 1, 1, (1 - z)/2), (4.70) 
and the asymptotic formula for F as n — > oo, we obtain: 



0u< n) 



L^^(^ +1 /2 + ^+l/2) . (4.71) 



dA s y^rm^x - y 

If we choose y c < x c , \y c \ < 1, we see that the critical value for x is x c = 1 and that 
the second term of the sum vanishes in the limit n — > oo. Defining the renormalized 
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length as I = na 2 / m , we obtain the expression (3.71). The computation of U2(p) 
and the fermionic loop operator v(p) can be done with the same method. We use 
several formulas for Gegenbauer polynomials: 



(xy)>C%\z) p - n - 6 z 



A3/2 ^ " V ^ ' ^ 



(%(z) = ^ n \ 2v ) , F ( n + 2u,-n,u+ 1 -, 1 -- Z -), (4.72) 

nV ; r(n + i)r(2i/) V 2 2 2 / 

and the asymptotic behavior of the Legendre polynomials. The final answers are 



(9 2 ± _ lu f dr 7 \ du 



at^ = -^ e " u {a-u +h )m- (474) 



where we have defined 



«(Z) 



and u(Z) = 2v / ttZu(Z). (4.75) 



We notice that these operators have an expansion at most quadratic in the fermionic 
couplings. This means that correlation functions involving more than two fermions 
vanish. In the next section we are going to derive the expressions for the scaling 
operators that follow from these loop operators. 
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5. SUPERSYMMETRIC MATRIX MODEL: 
ARBITRARY TOPOLOGIES 

In this chapter we will analyse the super-eigenvalue model in the double scaling 
limit and the connection to superintegrable hierarchies. One of the basic ingredi- 
ents to obtain the non-perturbative solution of the model are the super- Virasoro 
constraints in the double scaling limit. In order to calculate the non-universal parts 
of the loop operators we will compute explicit expressions for the scaling opera- 
tors on the sphere. We will find agreement with the scaling dimensions of (4m, 2) 
superconformal field theories coupled to 2D-supergravity. 

5.1. Continuum Super- Virasoro Constraints 

We will now see that, in the continuum limit, the bosonic and fermionic loop 
operators become a Z2-twisted scalar bosonic field and a Weyl- Major ana fermion 
in the Ramond sector. 

We redefine the purely bosonic part of the loop operator [18] according to: 

«o(p) = Mp) - v '(p) = - M (p)V(p- x )(p-y)- 

Using the relations (4.43), (4.44) and (4.47) we obtain in the scaling limit 

fc> 1 ^ + 2 ) OU 

This expression coincides with the Laplace transform of the loop operator. We 
will show that uq(z) can be identified with a free bosonic scalar field x(z) in two 
dimensions with antiperiodic boundary conditions 
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dx{z) =^a n+ iz- n -'i with x(e 2lvi z) = -x(z). (5.3) 
The Laplace transform of the bosonic loop can be decomposed into two pieces 

u Q (z) = t(z)+u (z). (5.4) 

Here t(z) denotes the non-universal part of the loop operator that is not dependent 
on the topology. The operator u = u$ + 112 can be expanded in terms of scaling 
operators (see [18] and references therein) in the following way: 



u(z) = k 2 J2 z k H°k) where (a k ) = ^. (5.5) 
fc>0 k 



Therefore, if we compute the Laurent expansion of the operator (4.62)* 

~¥ = ~Tz ~ g 4i L(l) \k + l + 2dt(U TdtT) ) Z 2 (5 - 6) 
for z — > 00, we obtain the following expression for the scaling operators {a k ) 



d(a k ) 1 (■ ) k T(k I - ! a 



dt k 2 k\V{\) \k+l 
The expansion of the macroscopic loop (4.73) 



+ 2d t (u k rd t r) ) . (5.7) 



* Now and in the following we redefine u(z), r and t in order to absorb — y c factors, and 
take effectively y c — 0. 
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(m = ^(-) k+1 y (^) - ^, where (a k ) = g, (5.8) 



determines the form of the scaling operators that are dual to the original couplings 
t k that differ from the couplings t k by a multiplicative constant: 



1 / n k+1 

d t(°k) = {j^ + idt^Tdtr) ) . (5.9) 



The relation between the coupling constants and is therefore 



h = H k+lI ^St k . (5.10) 



The non-universal part of uq(z) can be calculated using the planar string equation 
(3.85). After differentiation w.r.t. u we obtain 



d k+1 t 

__ = -t k+1 (k + l)\ + 0(u). (5.11) 



Now taking into account the relation between t k and t k we obtain the same result 
for the non-universal part of uq as in the purely bosonic model: 



t(z) = J2[k + ^)t k z k - 1 K (5.12) 
fc>0 ^ ' 

The bosonic-loop operator is therefore a bosonic scalar field with antiperiodic 
boundary conditions: 
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u(z) = (k + I) t k z k ^ + *-*- § te>- (5-13) 

k>0 ^ ' k>0 



The relations between the coupling constants of the model and the modes of the 
bosonic field are 



•9 , / 1\ t 



dt n ~ n_ 2 V 2 / K 



a n+ i = and a - n -- = [ n ~\ — ) — f° r n > 0. (5.14) 



We can handle in a similar way the fermionic-loop operator v(z) that will be 
identified with a Weyl-Majorana fermion in the Ramond sector 

^(z) = 5>„z- n -3. (5.15) 



First we redefine 



d(p)=v(p)-ap) = -^$- (5.16) 



With the recursion relations (4.50) we get, in the continuum limit, an equation 
analogous to (5.2): 



v(z)=a± ( -^ + 2 ]>>) fc -^^^ + ^ ) • (5.17) 



The scaling operators are calculated from the expansion of the loop operator and 
its Laplace transform 
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v(z) = K 2 a™ . 



n=0 



i °° l k 

i(v n ) and v(l) = Ka™J2(-) k+1 y(»k) 

fc=0 



(5.18) 



and are given by 



= ^^r# M V« and = -^ u V- (5 - 19) 

Note that the two-sided derivative in (5.19) guarantees Fermi statistics for the 
fermionic two-point functions 



9 ,~ ~ , 2 u d ,„ 

- {ViVn) = u _„« 



The relation between and is therefore 



(_)k+l Y(k + 

r k = 1 J,, U t 2 n- (5.20) 
fe 2A;! T(i) 

Using the planar equation r = X^x)^^' we ^ x ^ ne n on-universal part of the 
fermionic loop: 



^) = f 4 +E^ +l ' ( 5 - 21 ) 

Finally we obtain for the whole fermionic-loop operator a Weyl-Majorana fermion 
in the Ramond sector: 



97 



}(z) = ^z- 1 * + + K 2 J2 z ~ k ~ h ("k)- (5.22) 



fc>0 fc>0 



In terms of the mode expansion (5.15), we obtain the following representation of 
the coupling constants: 



*Po = 7T" + ^TT - ! i>n = kJt— and V'-n = — for n > 0. 

2/t aro OT n k 

(5.23) 

The identification of the zero mode guarantees = 1/2. The super-Virasoro con- 
straints in the continuum are therefore described by the super-energy-momentum 
tensor of a c — 1 superconformal field theory, 



T F (z) = ±dx(z)if>(z), (5.24) 

T(z) = i : dx(z)dx(z) : +)- : dx/>(z)i/;(z) : +-^. (5.25) 
With the mode expansion 



= ^ E ^""^n+i ^d T(*) = £ 2 L n , (5.26) 



we obtain in terms of the coupling constants: 



n+ 3 4« 2 ""' A V" ' 2 7 " K dr n+k+1 ' 2dt n '^ o ' K ^dt n+k+1 • ■' ^ dt k dr n . 

(5.27) 

The modes of T(z) are fixed by the anticommutation relation (1.48): 



ll.s 




dr n+ k+i 



d 




dtk^dtn^k 



d 2 



9 n / -i 

2 YV 2 



2/ dT k T n - k 



(5.28) 



After the double scaling limit we obtain that the superloop equations in the con- 
tinuum are equivalent to: 



with T(z,9) = T F (z) + 9T(z) determined by (5.24) and (5.25). 

The continuum limit of the superloop equations are described by a Z2-twisted 
massless scalar field and a Weyl-Majorana fermion in the Ramond sector. Notice 
that we only need to derive (5.24) since (5.25) follows from simple commutation 
relations. This makes the computation simpler than in the bosonic case, because 
Tp(z) has no normal-ordering ambiguities. 

Since we have explicit expressions for the scaling operators and the string 
equation for genus zero we can determine the critical exponents. For this purpose 
we have to consider the zero coupling regime, where t n = and rfe = 0. The 
expressions for the scaling operators follow from the free energy 



where D = d/dto. From the scaling of the free energy we obtain the string suscep- 
tibility T str as we saw in (2.43): 



Z _1 T(z, 9)Z = polynomial (z, 6), 



(5.29) 



k 2 D 2 F s =u- 2tD 2 t, 



(5.30) 
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dfF s ~ r Vstr . (5.31) 

We can read off the gravitational scaling dimension A of some scaling operator O 
from the relation (2.39) 

fi?(n°<)~ rr * tp+E<(A< " i) - ( 5 - 32 ) 

* i 

Using the form of the string equation for the m-multicritical point (3.85) we find 

dfF s ~u = t x l m . (5.33) 

We thus have the same family of susceptibilities T s t r = —1/m, m = 1, 2, 3, . . ., as 
the purely bosonic theory. For the bosonic operators we find 



dtM ~JT= U n d t U ~ t n/m+l/m-l = ^r str +n/m-l_ (5 34) 

ut n 

We obtain the same spectrum of dimensions as in the Hermitian one-matrix model: 



Tt 

K n = (5.35) 
m 

The scaling dimension of the fermionic operators can be determined from 

dti^n) ~ t k / m+n ' m - 2 = i-r st r+(A„ n -i)+(A, m -i)_ ( 5 3 ^ 

We obtain the result 
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These are the gravitational scaling dimensions of operators in the Neveu-Schwarz 
and Ramond sector of (4m, 2) N — 1 superconformal minimal models coupled to 
two-dimensional supergravity [41,42]. 



5.2. Iterative Solution of the Continuum Superloop Equations 

In the following we solve the continuum loop equations determined by (5.29), ex- 
panding the equations in the genus in a similar way as proposed by David for the 
bosonic theory [63]. In the continuum the loop equations get more transparent 
than in the discrete theory. It is interesting to consider also the planar case and 
see how the results for the loop operators and the string equation found in the 
previous sections appear in a simple way, as a solution of the continuum equations. 
For higher genera, i.e. for the torus, the double-torus, . . ., we obtain a systematic 
expansion determining all correlators beyond the planar limit. 

In the double scaling limit, the two superloop equations, that are equivalent 
to the continuum super- Virasoro constraints are: 





and 




v(z)dv(z) + k 2 (x W (z) + x FF (^) + A) = Polynomial (z), (5.39) 




where the two-loop operators are defined by: 
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y bf ( ) = y- z -k-i-2 d 2 F s 
x [ ) ^ dt k dn 

k,l>0 



dthdti ' 

k,l>0 



k,l>0 

The loop operators and the free energy have the following genus expansion 

u(z) = Uq{z) + KU\(z) + . . . = ^2 {^Q k \z) + KU\ k \z) + . . . j , 

fc>0 

v{z) = v (z) + KVl (z) + -.. = J2 ( v o k+1) ( z ) + ™f k+1 \z) + ...), 

k -° (5.41) 

X a (z) = xo(z) + kxHz) + ... with a = BF, BB, FF, 

F s = F + kFi + . . . . 

The sub-indices in our notation indicate the genus, while for the order in fermionic 
couplings we introduce upper indices. We now analyse the solution as a function 
of the genus. 

Planar solution 

The leading terms in the genus expansion of eqs. (5.38) and (5.39) are the planar 
loop equations 
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uq(z)vq(z) = polynomial (z), 



(5.42) 



uo(z) 2 — vo(z)d z vo(z) = polynomial (z). (5.43) 

We follow closely the steps of the discrete case with a one-cut ansatz for uq(z) and 
a similar expression for the fermionic-loop operator 



u (*) = M(z)V* + u + 7 and v (z) = r __ . (5.44) 



Expanding in powers of [z + u) 



JV(z) = J2 n k(z + u) k , M{z) = J2 m k(z + u)\ (5.45) 

fc>0 fc>0 

we see that A(z) is determined by demanding the r.h.s. of (5.42) to be polynomial 
in z 



Sl(!) = „,^-^^; (5.46) 



M(z) is determined from ~ 0(z 3//2 ), which gives 



dU ^ ~ ! ^-M(z) + (z + u)d u M(z)) - -^=^. (5.47) 



<9w ^/z + u \2 / 2y/z du 

The previous equation determines M(z) from the vanishing of the z~ l//2 contribu- 
tion 
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1 -M(z) + (z + u)d u M(z) = 1 -^. (5.48) 

Inserting (5.45) in this equation we obtain a relation between the coefficients rrik 
and the renormalized cosmological constant to: 



dt » a (-)fc r(i) d k+1 t 

mo = ^7 and mk = —¥(i^Y)^ for k>l, (5.49) 



and this fixes M(z) completely. The modes of N(z) are determined by demanding 
that vq(z) has for z — > oo the following expansion: 



v (z) = -^£L - V (z) = 0{z- 1 ' 2 ). (5.50) 
\JZ + u 

This implies the expansion dvo/du ~ 0(^ -1 / 2 ). Therefore we obtain the recursion 
relation 



°" k =-(k + -\ n k+1 for k > 1, (5.51) 



9m V 2 

that gives the value of all the n&'s it terms of no- Fermi statistics and compatibility 
between the bosonic and fermionic loops give the following relation between n\ and 
no: 

n\ = —Ad u nQ. (5.52) 

Therefore as a function of M(z) and uq the bosonic and fermionic loop operators 
have the following form: 
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m (z + u) 3 / 2 '" dt ^z + udt 

(5.53) 

With the ansatz no = Xm>o flnU n and (5.50) we obtain a relation to the fermionic 
coupling constants r n 



? (-mn + i) 

^n-r n - 2n!r( ^ r n . (5.54) 



Thus no is determined by r(w) in the spherical topology as follows: 

no = — t(u) with t(u) = r n u n . (5.55) 
The final result for the loop operators that follows from (5.53) is then 



tO+t 1 

uq(z) = M(z)\/z + u + and d t vo(z) = - dr. 

{z + u) 6 l z y/z + u 



This coincides with the results of sections 4.2 and 5.1. We obtain the planar string 
equation from the purely bosonic part of uq(z): 



M ° 0) -M(z)--^U = 0(z- 2 ). (5.56) 



\Jz + U \Jz + 



U 



The string equation (3.84) follows from the vanishing of the terms proportional to 
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fc>0 



Genus-one solution 



The two superloop equations obtained from (5.38) and (5.39) for genus one are^ 



2uqu\ - v\dvQ - vodvi + Xq B + Xq F + ^ = P°ly nomia l ( z ), (5.58) 



v\Uq + voui + Xo F = polynomial (z); (5.59) 

uq and vq are already determined since they are the solution on the sphere. The 
two-loop correlators are also determined by the results for genus zero and are 
obtained from (5.40) for Fg = Fq 



Y BF = j I 1 9 - 

2{z + uf/ 2 I ^T^tdt 



BB _ \ ( 1 _ M _ 1- f < (5 - 60) 

Xo "8\( Z + n) 2 z 2 J 2dt\(z + u) 



FF 1 / 1 1 



Xo 8V(^ + «) 2 *V 
We obtain for the fermionic-loop operator 



* In the following we will omit the z dependences. 
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v 1 = j2(z + u)- k - 1/2 V k+1/2 , 

k=0 



(5.61) 



with 



D 2 t\ „ DuDt 



The bosonic-loop operator can be determined in the same way. With the expansion 



4 

Ul (z) = J> + u)- k - 1/2 U k+1/2 , (5.62) 
fe=i 

we obtain for the coefficients 



D 2 u <-> 

U3/2 = Yon - 2tDV i/2 U 7/2 = 5rD(DuDr) 

12Du (5.63) 

U5/2 = — g- + 6^ 3/2J Dr - 2tDD 2 t U 9/2 = 7V 7/2 d t r. 

The value of V\j 2 cannot be determined by the requirement that the l.h.s. of 
(5.59) is a polynomial in z. It follows from the consistency between v\ and u\. 
This condition is given by the requirement of commuting derivatives: 



d 2 F 1 _ d 2 F 1 
dtdrk drkdt 



where F\ denotes the genus-one contribution to the free energy. 
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For the purely bosonic part of u\ we easily see that this is the same result as 
for the one-matrix model. The equation to be solved in this case is 



u (z)ui(z) + xo(z) + 



1 

8? 



polynomial (z), 



(5.65) 



Xo(z) = ^ ( , ,\ 2 + \ J 
8 \(z + uy z £ J 



(5.66) 



Equation (5.65) coincides with the order-zero equation in fermionic couplings com- 
ing from (5.58) just because the bosonic two-point correlators (5.60) and (5.66) 
coincide in the bosonic part. The reader may have noticed that the solution pre- 
sented for v\ is only linear in the fermionic couplings, as was the case for genus 
zero. This comes from the fact that the third order in fermionic couplings 



is solved by v \ = 0. This means that also on the torus we have a maximal 
coupling of two fermions. 

Genus-two solution 

Now we consider the situation in genus two to see if the heat capacity changes with 
respect to the one-matrix model. The two superloop equations to be solved are 

2uqU2 + u\ — V2&vq — vq8v2 — v\dvi + Xi* B + xf F — polynomial (z), (5.68) 



To calculate the heat capacity we are only interested in the purely bosonic piece 
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polynomial (z) 



(5.67) 



vqU2 + uqV2 + uiv\ + Xi — polynomial (z). 



(5.69) 



of the two-loop correlators that are given by (take (5.40) for F$ = F\): 



n Bb(°) 13 1 3 mi 1 (3m\ 5 m%\ 1 



We see that again these two-loop correlators coincide with the one-matrix model 
values. This immediately implies that, in genus two we have the same value in the 
purely bosonic part of the free energy as in the one-matrix model. The complete 
solution to order zero in fermionic couplings is 



(0) a b c d 



2 A 3 / 2 A 5 /2 ^ A 7 /2 ^ A 9 / 2 A n /2 

Q3mf 75m\rri2 145m 2 , 77mim^ 105m4 

~32m 7 + 16mg 128ml 32mJ~ + U8mJ 

Q3m\ 87mim2 105m3 

~ 32m[j ~ 32ml + 128m 4 

63m? 145m 2 , 203mi 105 
c=-— ^r + — -— r a~ 



32m 1 128m 4 128m 4 ""' 128m 3 / 

To compare this result with the one-matrix model we continue the calculation for 
the case of pure gravity (m^ = 0, k > 2). Our results reproduce up to genus two 
the values expected for the Painleve-I equation. Since 

1 2 

mo = and mi = -mf)D 2 u, (5-72) 

Du 3 

we get for the heat capacity the expansion 

( ^ o) = - 4 + 12* U7 J " 1~62 ( W J + • • • • (5 - 73) 
For pure gravity we take u = so that we obtain 
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1 ( r 1 k 2 49 k" \ ,. 
{aaao) = - I ^ t ----—— + ...y (574) 

This agrees with the first three terms of the solution to the Painleve-I equation 
appearing in pure gravity (3.123). 

5.3. non-perturbative solution of the discrete supersymmetric 
Model 

In the last section we have solved the superloop equations up to genus two. However 
the procedure is rather cumbersome and it would be desirable to find a more 
efficient method to obtain the solution for arbitrary genus. In this section we will 
analyse the discrete theory and find a powerful relation to the one-matrix model 
that will allow us to obtain the non-perturbative solution for the discrete theory 
[40]. 

We analyse the solution of the discrete superloop equations, or equivalently 
the discrete super- Virasoro constraints, through the representation (4.29) of the 
partition function. We first consider the case without fermionic couplings Ck+- = 
0. The partition function can be written as a function of the Pfaffian of the 
antisymmetric 2N x 2N matrix My-: 

M ij = \7' L = —L- for i^j, M ii = 

J Ai — Aj 

(5.75) 

and the usual Vandermonde A (Ay) = Yli<j 



r 2N ( 2N 2N \ 

Zs(9k,^ k+ l 2 = 0,27V) = / J]rfA n A(A ij )P/ 2iV (M ij )exp --^^^ • 

J n=l V S k i=l J 



(5.76) 
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Here the product of the Vandermonde and the Pfaffian is totally symmetric under 
exchange of two eigenvalues. We would like to show that the relation of the su- 
persymmetric partition function at zero fermionic couplings (5.76) to the bosonic 
partition function of the one-matrix model (3.1) is given by: 

Zs(9k,C k+ i = 0,27V) = ±j^Z 2 B (9k,N). (5.77) 

To prove this equality we show that the following equivalence for the measures of 
both models holds 

2N 

Pf N {M l3 ) Y[(\i - Xj) = W J2 S J ' J A 2 (/) A 2 ( J). (5.78) 

i<j I, J 

Here we have introduced the notation / = (A^, . . . , \i N ), i\ < ... < ijv (same 
for J), and S I,J means symmetric permutations between the elements of / and 
J. The equality (5.77) is then a consequence of (5.78) since the partition function 
(5.76) with 2N eigenvalues will factorize into ( 2 ^) products of two independent 
one-matrix models with N eigenvalues. 

Obviously (5.77) holds for the case of two eigenvalues, so that to illustrate the 
situation we start with the first non-trivial case of four eigenvalues N = 2. We will 
consider the left- and r.h.s. of (5.78) as polynomials -P(A) and Q(X), respectively, 
in A = A4: 

P(A) = (A - A 2 )(A - A 3 )Ai 2 Ai 3 - (A - Ai)(A - A 3 )Ai 2 A23 + (A - Ai)(A - A 2 )Ai 3 A 2 3, 

(5.79) 

Q(A) = i ((A - A 3 ) 2 A? 2 + (A - A 2 ) 2 A 2 3 + (A - A!) 2 A 2 3 ) . (5.80) 
Because both polynomials agree for A = A^, % — 1, 2, 3, and since they are of degree 
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two in A, it follows that Q(X) = -P(A) for all A. 

We show that (5.78) is correct by induction over N, i.e. assuming that (5.78) 
holds for 2N eigenvalues we prove the case of 2N + 2 eigenvalues. For Af > 2 we 
set A = A2A/+2; Q(X) and P(A) are then polynomials of degree 2N — 2 in A. To 
show that both polynomials agree, we prove that they agree in 2N + 1 points. It 
is enough to show that they coincide for A = \2N+i since the symmetry under 
the interchange of \2N+i and for i = 1, . . . , 2N guarantees -P(Aj) = Q(\i) for 
i — 1, . . . , 2N. Using the induction hypothesis we see that Q(X 2 n + i) agrees with 
P(X 2N+1 ): 

Q(A 2 iv + i) = E S J ' J A 2 (/)A 2 ( J) Y[ (A, - X 2N+1 ) 2 (X j - X 2N+1 ) 2 

I, J h<...<iN 
ji<...<j N 

2N 

= Pf 2N (M tj )l[(X i -X J ) J] (Ai-A 2iV+1 ) 2 J] (X J -X 2N+1 ) 2 = P(X 2N+1 ). 

i<j h<...<iN ji<—<jN 

(5.81) 

This completes the proof that Q(X) = P(X) in the case of 2iV + 2 eigenvalues. We 
conclude that the bosonic part of the supersymmetric 2N x 2N model is propor- 
tional to the square of the corresponding N x N bosonic one- matrix model (5.77). 
With this result we can already find the form of the discrete string equation for 
arbitrary genus. 

Writing Z = e N2ps , eq. (5.77) implies that the fermionic independent piece of 
the supersymmetric free energy is related to the free energy of the one-matrix 
model J-g (up to an irrelevant additive constant): 

F$» = f . (5.82) 
Since the dependence of the free energy on go is trivial, dF/dgo = — 1/A, (5.82) 
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implies a rescaling of the cosmological constant A5 = 2Ag. This relation can also 
be obtained, if we evaluate the bosonic piece of the Lq constraint (4.25): 



A s = 2j2kg k d As -A|y^ • (5.83) 
fc=i ^ 9k ' 

Up to a factor of two, the cosmological constant satisfies the same string equation 
as in the one-matrix model (3.39). 

Including the fermionic couplings, the following equality for the second order 
in fermionic couplings holds: 



d 2 Z s (2N) 1 (2N\ ( d 2 Z B {N) dZ B {N) dZ B {N) \ 

dZ k+ .dZ n+ .-2»AN){d gk+1 d gn ZB{N} d 9k+1 dg n ) {k ~ n) - 

(5.84) 

This formula is purely bosonic in the sense that after derivation we set the fermionic 
couplings to zero. The simplest example is N = 1 



f Z f ] = I dX.dX^M - A 2 )(A?A* - A{AS)e*> (Af+AS) 

= j d\id\ 2 (\ k+1 \% + a£ +1 A? - A? +1 A£ - A2 +1 Af) e 5>( A * +A 



\dg k+1 dg n dg k+1 dg n J 



(5.85) 



The general proof can be done with similar methods as the previous ones. We have 
to take into account the identity: 



113 



2N r 2N 

a,/3=l J i=l 



- 5 J ' J A 2 (J) A 2 ( J) (^K - X t lx D "(*<-")■ (5.86) 



I,J a,/3=l 



This relation can be proved by comparing a general coefficient Xi U \j from the left 
and r.h.s. of the above equation. We have to use the property (5.78). Equations 
(5.82) and (5.84) imply that the free energy, until the second order in fermionic 
couplings, is given by: 



Since the solution to the continuum super- Virasoro constraints will be analysed 
in a moment, it will become clear that (5.87) is the complete expansion of the free 
energy and no dependence on more than two fermionic couplings appears. This 
characteristic has been confirmed by McArthur [71] using methods similar to the 
previous ones. Equation (5.87) gives a representation of the model (4.29) in terms 
of the one-matrix model since the fermionic variables 9{ have been integrated out. 

5.4. Solution to the Super- Virasoro Constraints in the Double 
Scaling Limit 

In section 4.1 we saw that the super- Virasoro constraints are described by the 
super-energy-momentum tensor of a c = 1 superconformal field theory. The modes 
of the supercurrent Tp = \ *Yl<G n+ i_z~ n ~ 2 are described in terms of the coupling 
constants of the model for n > —1: 




(5.87) 
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+ ' 4 « 2 ""' " ' ^ V" ' 2 7 ""^m+fc+i ' 2 dtn ' t^ ' W< ^n+fc+i ' " ^ dt kdr n - k 

(5.88) 

The constraints imposed on the partition function Zs(t k , T k ) = exp^s) are: 



G n+ iZ 5 (t fc ,r fc ) = for n > -1. (5.89) 

We now present the solution to these constraints in the double scaling limit. The 
constraints determine that in this limit the model (4.29) is described by the free 
energy 



k 2 DF s = D~ l u - 2tDt, (5.90) 

where D = d/dto and t is the renormalized cosmological constant satisfying the 
string equation 



t = - ^2(2k + l)t k R k [u]. (5.91) 
k>\ 

The non-perturbative form of the fermionic scaling variable r = dF s /dro — ro/2 is 

r = -J2r k RkH (5.92) 

fc>0 

As usual, the bosonic piece of the two-point function of the puncture operator u 
satisfies the KdV-flow equation 
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— = DR n+1 [u], 



(5.93) 



where Rk[u] are the Gel'fand-Dikii polynomials. 

To show that eqs. (5.90), (5.91) and (5.92) are the solution to the super- 



constraints are satisfied, since the algebra (1-48) guarantees for all the remaining 
constraints. That this is indeed the case can be shown by simple calculations, where 
we have used well-known relations of the one-matrix model, like the recurrence 
relations of the Gel'fand-Dikii polynomials. 

This form of the free energy, the string equation (5.91) and the form of r (5.92) 
coincide with the solution presented in sections 4.2 and 5.2 for the first few orders 
of the genus expansion. The free energy (5.90) can be written as a function of the 
free energy of the one-matrix model; this is very similar to the relation (5.87) of the 
discrete model and has an expansion at most quadratic in the fermionic couplings: 



The form of the free energy and the string equation obtained imply that the 
bosonic piece of the model is equivalent to an ordinary one-matrix model in the 
double scaling limit, where the constraints act on the square root of the partition 
function L n y/Zj$ = for n > — 1 [18]. 

The fermionic dependence of the free energy determines that correlation func- 
tions involving more than two fermionic operators vanish. These properties were 
known to hold for the first orders of the genus expansion, as shown in sections 4.2 
and 5.2, and are in agreement with the previous analysis of the discrete model. 



Virasoro conditions (5.89), it is enough to show that the G_i/2 an d the G3/2 




(5.94) 



k>i 
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The correlation functions in the double scaling limit follow from (5.90), (5.91) 
and (5.92). As a function of the two-point function of the puncture operator 
U = u — 2tD 2 t, they are given by: 



WT k = 2D H"lBr) • f = DR M M - 2 A (tD\) . (5.95) 



The scaling dimensions of the scaling operators (5.95) and the string suscep- 
tibility obtained from (5.91) are in agreement with the results of (4m, 2)-minimal 
superconformal models coupled to 2D-supergravity [42,41]. From (5.95) we obtain 
for the first non-trivial bosonic flow a set of two equations: 



— = D(3u 2 + k 2 D 2 u), — = 6uDt + k 2 D\ (5.96) 



where the first one is the ordinary KdV equation. These equations are invariant 
under the global supersymmetric transformations 



5u = eDr and 5r = eu, (5.97) 

where e is a constant anticommuting parameter. This property is known to hold for 
the supersymmetric extensions of KdV of Manin-Radul [29] or Mulase- Rabin [72] . 
We will introduce the supersymmetric extensions of the general KP hierarchy in the 
next section and restrict ourselves for the moment to describing some properties of 
the correlation functions. We can write the form of our generalized KdV equation 
in terms of U the variable that is, perhaps, more convenient: 
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= D(3U 2 + k 2 D 2 U - UDUtDt + Qk 2 DtD 3 t), (5.98) 



— = WDt - 12tDtD 2 t + k 2 D 3 t. (5.99) 

Equations (5.98) and (5.99) can independently be derived from the first orders of 
the genus expansion, with the method proposed in ref. [37]. The bosonic flows 
satisfy the following recursion relations that can be obtained from the recursion 
relations of the i?fc[u]'s: 



du/dt k+x \ / k 2 D 2 + 2u + 2DuD~ 1 
dr/dt k+1 J ~ \ 2DtD- x + 2D- 1 Dt 

The fermionic flows satisfy the relations: 



\ (du/dt k \ 

k 2 D 2 + 2u + 2D- x uD J \dr/dt k J 

(5.100) 



-^=0 and = (k 2 D 2 + 2u + 2D- l uD)^. (5.101) 

or k dr k+1 dr k 

The aim of the next section is to introduce the supersymmetric extensions of 
the KP hierarchy that are known in the literature, and to see whether it is possible 
to find agreement with our correlation functions. 

5.5. Super-Integrable Hierarchies and 2D Supergravity 

In the last section we obtained the form of the correlation functions of scaling oper- 
ators in the double scaling limit. To gain insight into the geometric interpretation 
of the solution obtained it would be interesting to see whether we are able to iden- 
tify the hierarchy of non-linear differential equations obtained from our model with 
one of the supersymmetric extensions of the KdV (or more general KP) hierarchies 
known in the literature. 
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The program of supersymmetrization was initiated by Manin and Radul [29]. 
They introduced a supersymmetric generalization of the KP system in the Lax 
formalism. The Manin-Radul hierarchy [29] is a set of flow equations deforming 
the pseudo-differential operator of the form: 



Q = 5 + Q + Qi5- 1 + Q 2 5- 2 + ... (5.102) 

where 



*=^ + %, ("OS) 

and the coefficients Ui(x,£,T) are formal power series in the even variables x : T2n 
and the odd variables £,T2 n _i, with n — 1,2, . . .. We will use the notation T = 
(Ti,T2, ...). The expression for can be obtained using 5 _1 = Sd^ 1 , where we 
have taken into account that 5~ 2 = d^ 1 . 

The condition 5Qo + 2Qi = is necessary and sufficient for the existence of a 
pseudodifferential operator: 



S = 1 + SkT 1 + S 2 5~ 2 + . . . with S~ 1 QS = 5. (5.104) 
In terms of the operator Q, the Manin-Radul hierarchy reads: 



= [Q2» Q] = _ [Q 2» g] (5 1Q5) 

oT 2n 

= -[Q 2 - n -\Q] + f>-x^ ^— , (5.106) 

where Q n = Q\ + Q- is the decomposition into non-negative and negative powers 
of 5. We can also describe the hierarchy in terms of the operator S: 
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Ol 2n 



(5.107) 



S. (5.108) 

The term involving the infinite sum was necessary to make the Manin-Radul system 
completely integrable, but this term also made it quite difficult to understand its 
geometric meaning. 

To understand whether there is a relation between our correlation functions 
and this hierarchy, it is important to know the r-function. This super- r-function 
was computed by Ueno and Yamada [73] for the Manin-Radul hierarchy in terms 
of superdeterminants. The partition function of the matrix model can then be 
identified with the r-function of the hierarchy, so that we have a unique way of 
identifying the variables from the matrix model with those of the superhierarchy. 
Di-Francesco et al. [28] have analysed the role that this super-r-function plays in 
the description of N — 1 unitary SCFTs coupled to 2D-supergravity. Using the 
Lax-pair formalism they identified the string equation and the relevant primary 
fields. However, they came to the conclusion that odd flows are inconsistent with 
the string equation, so that it does not seem so clear that this hierarchy plays a 
role in 2D-supergravity. 

A different supersymmetric generalization of KP, that precisely differs from the 
previous one in the odd flows, is the Mulase- Rabin hierarchy [72]. It is formulated 
in terms of the operator S rather then in terms of a Lax pair. It is a more natural 
supersymmetric generalization of ordinary KP, from the geometric point if view. 
The flows are given by: 



dS 



0T 2 



S 5 



2n-l 



:n-l 



oo 

k=l 



2k- 



^2n+2fe- 



OS 



■[ss^s-^s 



(5.109) 
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-^±— = -[sdcd^s-^s = -[s^- 1 - tf^s-^s. (5.110) 

<J±2n-\ 

If we compare these equations with the Manin-Radul hierarchy in (5.107) and 
(5.108), we realize that the even flows are the same while the odd flows are com- 
pletely different. Both hierarchies share the characteristic that if the fermionic 
variables are set to zero we obtain the standard KP hierarchy. This is a common 
feature of all the supersymmetric extensions of KdV (KP). 

Since the above formulation is rather abstract, we would like first to compute 
explicitly the simplest even and odd flows of the most general two reduction. The 
sKdV equations, i.e. the supersymmetric generalization of KdV, are obtained by 
considering the flows obtained from the operator [29] : 



L = 5 4 + v\5 + vo, (5.111) 

where v\ is fermionic and vq bosonic and L = Q A . This is the most general two 
reduction of both supersymmetric hierarchies, since the definition of S in terms of 
Q (5.104) determines v 3 = i>2 = 0. To calculate the flows we have to take into 
account the generalized Leibniz rule: 

5{ab) = 5(a)b+(-fa5(b), 

3=0 

Here a = if a is a bosonic function and a — 1 if a is fermionic. From the definition 
of S in terms of L, 



(5.112) 
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LS = S5 4 , (5.113) 
we can obtain a representation of S in terms of v i and vq: 



S 2 S 1 = 

2 ' 



A2c VlSl V ° rrttA\ 

S 2 = — —, (5.114) 



5 2 S 3 = ^ - ^M-Si + + uqS'i). 



4 2 

For arbitrary we have the recursion relation 



(-) k+1 1 

S 2 S k+1 = y -^—v 1 S k - -{5 A + v 1 5 + v )S k _ 1 for fc = 2,3,4,... (5.115) 

These relations are valid for both the Manin-Radul and the Mulase-Rabin hierar- 
chies, since we only need, to obtain them, the definition of S in terms of L and not 
the special form of the flows. 

The first bosonic flow for these hierarchies is: 



dv dv dv\ dv\ , 

w 2 = ^ and wr^- (5 - 116) 

These are the equations for a chiral wave moving to the left. A similar equation 
appears also in the ordinary KdV hierarchy. The flow with respect to T4 is trivial, 
taking into account the special form of the operator that we are considering (5.111). 
The next even flow leads us to the generalization of the KdV equation to the 
supersymmetric case 
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8vq d ( 1 8 2 vq 3 2 , 3 



+ 7^0 + 7^0 , (5.H7) 



dT 6 dx V 4 dx 2 4 u 4 



dvi d (\ d 2 vi 3 _. 3 . _ , , 

+ -vidvi + -v\vo . (5.118) 



<9T 6 &r V4 dx 2 4 2 

Since the Manin-Radul and the Mulase-Rabin hierarchies agree in the bosonic 
flows, these equations are valid for both systems. Setting vi — we get the 
ordinary KdV equation. Setting vq = and vi — wi + £wo, where ^ is a purely 
fermionic parameter, (5.117) takes the form: 



dwo d ( 1 8 2 wq 3 2 3 dw\ 



dT 6 dx V 4 dx 2 4 u 4 dx 



+ 'A ~ 7^1^T . (5-119) 



dwi d (\ d 2 wx 

+ -wqWi . (5.120) 



dT 6 dx\4 dx 2 4 

Setting wi = we again get the KdV equation. This means that there are two 
possible operators that will lead us to the KdV hierarchy: 



Li=5 4 + £w 5 and L 2 = 5 4 + v . (5.121) 

In general, it is not consistent to consider the odd flows of these operators. The 
coefficients v\ and wi are uniquely determined by the r-function, so that we have 
no freedom to set them to zero and still keep odd flows consistently in the theory. 

In the following we will write down the expressions for the simplest fermionic 
flows of the operator (5.111). First we consider the Mulase-Rabin case. The form 
of the first fermionic flow is then given by 
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dSi dSi , dS 2 Q _ Q dS 2 ,_ 100 , 

7^ = ^ and — = -S 3 + S X S 2 + — . (5.122) 

The second fermionic flow satisfies 



H = 5 (S 2 S 1 + 5S 2 + S 3 - S1S2) - (^i) 2 . (5.123) 

The fermionic flows are local in terms of the variables Si, but not in terms of v\ 
and vq. For L\ we can write the second fermionic flow in the following way: 



^r = \ [(S-^wo) Stw + 6 3 tw Q } . (5.124) 

This equation can be expanded in the number of fermionic couplings. The purely 
bosonic part of the resulting equation leads us to 5 2 wq — wq5~ 2 wq = 0, which 
clearly is inconsistent. Thus we cannot allow the operator L\ as a reduction of 
Mulase- Rabin and take the odd flows into account at the same time. The operator 
L = <5 4 + v \S was considered in ref. [74], with a similar conclusion. 

The first fermionic flows of the Manin-Radul hierarchy are: 



or gi^-g^ ( , 125) 



^ = 5 (5 2 S! + 5S 2 + S 3 - SiS 2 - Si5Si) T%.A. (5.126) 

Ol 3 ^ Ol 2k+2 

The expressions for dSi/dT^ are easy to write down because the bosonic flows 
for L\ have the same expressions as for ordinary KdV. More general flows can be 
obtained (in principle) using recursion relations [75,74]. 



124 



As we said previously, to see if one of these hierarchies describes our correlation 
functions in the double scaling limit we have to find the identification between the 
variables of the supersymmetric model with those of the hierarchy. The easiest way 
will of course be the identification of the r-function with the partition function. The 
first even and odd flows of the Manin-Radul hierarchy in terms of the r-function 
take the form [73] 

fl(logr) = (Ji(logr) = Si, 
5 2 (logr) = -^ 1 , 



<5 3 (logr) = - Q(35 + 5i)S 2 + d x S! - SrfS^j , 



5 4 (logr) = 25S 3 + SdxSi - (5S!) 2 - 26(S 1 S 2 ). 
Here we have used the notation 



hi 



k=1 --2l+2k-2 



(5.127) 



(5.128) 



It is easy to show that the following identities hold: 

SiiSr) = -6(S!) 

Si(S 2j ) = -5(S 2j ) - 2S!S 2j + 2S 2j -i (5.129) 
h(S 2 j-i) = S(S 2 j-i) - 2S\S 2 j-i. 
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Eqs. (5.127) give a relation between the coefficients vq and v\ of the operator 
L (5.111) and the r-function. The first equation implies the condition v\ = 
— 25 3 (logr). Requiring compatibility among both KdV equations fixes t>o = 2<5 4 (logr).| 
This means that the operator we are going to consider in the comparison with the 
Manin-Radul hierarchy is 

L = 5 4 - 25 3 (logr)5 + 25 4 (logr). (5.130) 

If we now postulate that the partition function is the r-function of the KdV hier- 
archy, this would mean that (up to multiplicative factors coming from the identi- 
fication of the coupling constants of the matrix model) the two-point function of 
the puncture operator 

U = u-2tD 2 t (5.131) 

satisfies the KdV-flow equations. 

To see if we can find agreement with the matrix model, we can check, for 
example, if the KdV equations (5.117) and (5.118) are satisfied. We already know 
how to interpret vq and v\ in terms of correlation functions. If we now check the 
KdV equations we arrive at the conclusion that the fermionic dependence in these 
equations is inconsistent!. 

We have no other choice for the identification of the variables between the two 
theories. We conclude that the Manin-Radul hierarchy does not reproduce the 
flows of the super-eigenvalue model. 

This inconsistency will appear also in a comparison with the Mulase-Rabin 
hierarchy, since we are considering an even flow. However in this case we have 
the difficulty that there does not exist a formulation in terms of a r-function yet 
[76], so that in principle the connection between vq and v\ and the matrix-model 
variables could be different. 
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The above identification is the simplest generalization we could think of, since 
it is a formulation in terms of the two-point function of the puncture operator. 

Next, we would like to remark that there exists another extension of the KdV 
hierarchy that has been proposed by Kupershmidt [77,78]. In this hierarchy the 
KdV equation takes the following form: 

d 

u t = — (3u 2 - u xx + 3if(p x ) (5.132) 



<f t = 3u x <p + 6u<f x - 4if xxx , (5.133) 

where u is a commuting and if is the anticommuting variable. We notice that 
while the first equation coincides with the first one of Manin-Radul (after a trivial 
rescaling), the second one is different. In analogy to the bosonic KdV hierarchy, 
this one has a formulation in terms of two Hamiltonian structures and a non-trivial 
Lax representation. However this system of differential equations is less interesting 
for the comparison with the correlation functions of the super-eigenvalue model, 
since the odd flows are not present in this formulation. 

We therefore take a more radical point of view and claim that the hierarchy 
of super-differential equations presented is a new supersymmetric extension of the 
KdV hierarchy. The properties of this hierarchy have been studied by Figueroa- 
O'Farrill and Stanciu [79] and by McArthur [80]. In ref. [81] a discrete analogue 
to this integrable hierarchy, that is related to the Toda lattice, is found. 

It would be interesting to see whether it is possible to find a connection between 
this super-hierarchy and topological supergravity in the same sense as it holds 
for the one-matrix model [18]. Furthermore it would be interesting to find the 
connection to the Landau-Ginzburg approach presented by Di Francesco et al. [28] 
and to see whether the new hierarchy explains the discrepancies found. 
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6. CONCLUSIONS 



The (discrete) Virasoro constraints, which are closely related to the integrability 
of the one-matrix model, can be generalized to the N — 1 supersymmetric case in 
terms of the modes of the energy-momentum tensor of a c = 1 free massless super- 
field. Imposing the super- Virasoro constraints on the partition function, Alvarez- 
Gaume et al. [30] have constructed a discrete super-eigenvalue model formulated in 
terms of N even and N odd eigenvalues (Aj, 9{) and bosonic and fermionic coupling 
constants (gk,£k+i/2) respectively. This model can be solved by evaluating the su- 
perloop equations, which nicely generalize the loop equations found by Kazakov 
for the one- matrix model [14]. 

As a starting point, we have solved this model for generic potentials and genus 
zero by evaluating the planar superloop equations. The solution led to a set of 
multicritical points described by a planar string equation that has the same form 
as in the one-matrix model. These multicritical points are labeled by an integer 
m = 1,2,..., where the string susceptibility takes the value T str = —1/m. We 
have obtained explicit expressions for the macroscopic superloop operators and the 
microscopic scaling operators, (a n ) and (v n ), corresponding to the operators in the 
Neveu-Schwarz sector and Ramond sector of the SCFT coupled to 2D-supergravity. 
The values of the gravitational scaling dimensions and the string susceptibility 
correspond to those of (4m, 2)-minimal SCFTs coupled to 2D-supergravity [41,42]. 
The correlation functions of operators in the Neveu-Schwarz sector agree with those 
computed in the super-Liouville approach for spherical topologies [8,31,32]. 

Considering arbitrary potentials was important in order to derive the Virasoro 
constraints in the double scaling limit. These constraints are formulated in terms 
of a c = 1 theory with a Z2-twisted scalar field and a Weyl- Major ana fermion in 
the Ramond sector. The superloop equations in the continuum limit are equivalent 
to the double scaled super- Virasoro constraints. 

We have solved these superloop equations to all orders in the genus expan- 
sion. Up to genus two we have used the expansion of these equations in the string 
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coupling constant to obtain the form of the microscopic scaling operators. For ar- 
bitrary genera we have been able to express the free energy of the super-eigenvalue 
model in terms of Gel'fand-Dikii polynomials. Using their recurrence relation it is 
possible to show that the super- Virasoro constraints are satisfied. It turned out, 
that the bosonic part of the model (fermionic couplings set to zero) reproduces 
the solution of the one-matrix model. We have also been able to write the super- 
eigenvalue model as a matrix model in terms of purely bosonic hermitian matrices. 
The non-perturbative expression for the free energy of the super symmetric theory 
had a simple connection to the free energy of the discrete one-matrix model. This 
relation is important in order to find the geometric interpretation of this model 
and to obtain the Feynman rules for the discrete theory. 

The non-perturbative solution of the continuum theory showed that the inte- 
grable structure that appears behind this model is a new supersymmetric general- 
ization of the KdV hierarchy. 
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